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A new model for space and matter is obtained by joining every pair of point charges in the 
observable universe by an ethereal string. Positive gravitational potential energy in each string gives 
an attractive gravitational force due to the action of an energy conservation constraint. Newton’s 
laws of motion are derived and inertia is explained in accordance with Mach’s principle. The 
Machian string model gives a surprisingly simple way to understand the expansion history of the 
Universe. The decelerating expansion in the radiation era and the matter era is explained without 
using General Relativity and the transition from deceleration to acceleration is explained without 
the need to introduce a separate ‘dark energy’ component. The interaction between Machian strings 
gives a physical model for modified Newtonian dynamics (MOND) and is therefore an alternative 
to ‘dark matter’. 

“ I consider it entirely possible that physics cannot be based upon the field concept, that is on continuous 
structures. Then nothing will remain of my whole castle in the air, including the theory of gravitation, but 
also nothing of the rest of contemporary physics.” (Albert Einstein) 

PACS numbers: 


I. MOTIVATION 

The purpose of the present paper is to suggest a new 
model for the structure of elementary particles and a new 
model for space and to offer a new approach to the prob¬ 
lems of ‘dark matter’ and ‘dark energy’. 

One of the basic properties of elementary particles is 
that they have spin one-half. Quantum mechanically, a 
spin one-half particle ends up in a different state when ro¬ 
tated through an angle of 27 t and only returns to its orig¬ 
inal state when rotated through an angle of 47 t. There is 
no analogue of such behaviour for a classical point par¬ 
ticle but a mechanical model for spin one-half can be 
constructed when the particle is connected to its sur¬ 
roundings, for example by someone holding onto it or 
by strings connecting it to a rigid external frame mm- 
Suppose that the strings are initially not entangled. A 
rotation of the particle through an angle of 4-7T results in 
a configuration that is equivalent to the initial configu¬ 
ration, in the sense that the strings can be disentangled, 
but a rotation through an angle of 27 t results in a config¬ 
uration that is not equivalent to the initial configuration. 

In searching for a new model for space it is natural to 
reconsider Mach’s principle, which was never fully incor¬ 
porated into General Relativity I3H5]- Instead of defining 
inertial frames in terms of a local absolute space, Mach 
argued that the inertial forces experienced by an acceler¬ 
ated particle are due to some sort of interaction with the 
distant matter in the Universe. Perhaps a physical model 
for space can be constructed by defining space itself in 
terms of distant matter [6]. 

In the model of Machian space quanta, the mechanical 
model for spin one-half and Mach’s principle are incorpo¬ 
rated by postulating that an elementary particle consists 
of a set of ethereal strings connecting a charged centre 
to the centre of every other elementary particle in the 
observable universe. The set of strings associated with a 
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given elementary particle can be imagined as a quantum 
of space with the size of the observable universe and will 
be referred to as a Machian space quantum. The strings 
contain energy, and hence mass, but are not charged. The 
strings associated with the space quanta of two nearby 
particles of mass mi and m 2 are illustrated in Figure [I] 
The strings connecting two such masses will be referred 
to as direct strings and strings connected to distant mat¬ 
ter will be referred to as Machian strings. 



FIG. 1: Schematic diagram of the direct strings and the 
Machian strings in the Machian space quanta of two masses 
mi and m 2 ■ The direct strings join the centres of the space 
quanta associated with the two masses and the Machian 
strings connect the centres with the centres of all the other 
space quanta in the observable universe. 

The idea that an elementary particle could have the 
size of the observable universe appears, at first sight, to 
be completely contrary to experience. Indeed, high en¬ 
ergy scattering experiments have failed to find any spatial 
extent at all in the charge distribution of either quarks or 
electrons. However, experiments cannot probe an entire 
space quantum because the only part of a space quantum 
that can interact directly with another space quantum 
is the centre. The only conclusion that may therefore 
be drawn from scattering experiments is that the charge 
distribution at the centre of a space quantum is pointlike. 

The gravitational interaction between the two masses 
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mi and m 2 illustrated in Figure [l] is defined by speci¬ 
fying the total energy, Gmim^/r, in the direct strings 
joining them, where r is the distance between the centres 
of the two masses. To a first approximation, the gravita¬ 
tional masses m\ and m 2 are simply the total masses in 
the associated space quanta and it makes no difference 
at all whether all the mass of a space quantum is at the 
centre or whether the mass is distributed throughout the 
strings. For the further development of the theory, how¬ 
ever, the fact that the mass is distributed throughout the 
strings is very significant. The fundamental variables in 
the problem are the strings themselves, rather than posi¬ 
tions and velocities of the particle centres. The energies 
in the strings may be calculated explicitly, including the 
dependence on the velocities of the particles, and it is 
then possible to derive the properties of particle inertia 
and kinetic energy from first principles instead of putting 
them in by hand. When the contribution to the energies 
of the strings from the expansion of the Universe is taken 
into account, the string model is found to give an accel¬ 
erating expansion of the Universe without the need to 
introduce any mysterious ‘dark energy’. 

It might be thought, at first sight, that Machian strings 
give a very simple solution to the problem of ‘dark mat¬ 
ter’. It is well known that flat galaxy rotation curves can 
be accounted for by assuming that galaxies lie within 
a dark matter halo whose density at a distance r from 
the centre of the galaxy is proportional to 1/r 2 [2- For 
strings with uniform mass per unit length, the mass in 
the strings enclosed within a spherical shell of thickness 
dr at any radius is proportional to dr and it follows that 
the density in the strings is proportional to 1/r 2 , which 
appears to be exactly what is required. This observa¬ 
tion was, in fact, the original motivation for the Machian 
string model. Unfortunately the idea doesn’t work be¬ 
cause the only gravitational interaction between space 
quanta is the interaction due to the energies in their 
common direct strings. Unlike a conventional mass dis¬ 
tribution, the mass elements within a Machian string are 
invisible to other space quanta and do not, therefore, act 
as additional point sources for the gravitational field. In 
any case, such a model would give an additional gravi¬ 
tational acceleration proportional to M/r, whereas the 
additional acceleration required to account for galaxy ro¬ 
tation curves is proportional to y/~M/r. A different model 
for ‘dark matter’, based on a contact interaction between 
the strings of different space quanta, is introduced in Sec¬ 
tion 0 

Whereas the quantisation of matter is automatically 
present at the classical level a photon is, by definition, a 
quantum mechanical object. The strange quantum me¬ 
chanical behaviour of the photon was described by Ed¬ 
dington in 1928 as follows. 

The pursuit of the quantum leads to many sur¬ 
prises; but probably none is more outrageous to 
our preconceptions than the regathering of light 
and other radiant energy into h-units, when all 
the classical pictures show it to be dispersing 
more and more. Consider the light waves which 
are emitted as a result of a single emission on the 
star Sirius. These bear away a certain amount 
of energy endowed with a certain period, and the 


product is h. The period is carried by the waves 
without change, but the energy spreads out in 
an ever-widening circle. Eight years and nine 
months after the emission the wave-front is due to 
reach the earth. A few minutes before the arrival 
some person takes it into his head to go out and 
admire the glories of the heavens and - in short 
- to stick his eye in the way. The light waves 
when they started out could have had no notion 
of what they were going to hit; for all they knew 
they would be bound on a journey through end¬ 
less space, as most of their colleagues were. Their 
energy would seem to be dissipated beyond recov¬ 
ery over a sphere of 50 billion miles’ radius. And 
yet if that energy is ever to enter matter again, if 
it is to work those chemical changes in the retina 
which give rise to the sensation of light, it must 
enter as a single quantum of action h [8j. 

How can a photon possibly manage to keep itself all to¬ 
gether as a single tiny quantum when it is so incredibly 
spread out across the Universe ? The answer, perhaps, is 
that a photon is also a Machian space quantum. A pho¬ 
ton space quantum still has a centre, defined as the point 
of intersection of all the photon strings, but the centre 
has no mass and no charge. 

Machian space quanta may also help to understand 
the invariance of the speed of light. Special Relativity 
is based on the postulate that the speed of a light wave, 
c, is independent of the motion of the light source rel¬ 
ative to the observer. Is there any way to understand 
the invariance of c with some sort of physical model or 
must the postulate simply be accepted ? In the old ether 
model, the speed of light was assumed to be equal to c 
in the rest frame of an ‘ether’ that pervades all space. 
The invariance of the speed of light waves with respect 
to motion of the source is then obvious, just as the speed 
of sound waves in air and water waves in water are also 
independent of the motion of the source. However, the 
invariance of c with respect to motion of the observer 
was much harder to explain and led to the conclusion 
that motion through the ether somehow causes the time 
dilation of clocks and the length contraction of metre 
sticks. The Machian string model offers an alternative 
point of view because the strings of the observer space 
quantum define an ether that is always at rest relative 
to the observer. Suppose the speed c is a characteristic 
of the interaction between a photon space quantum and 
the space quantum of an observer. The invariance of c 
with respect to motion of the observer then follows auto¬ 
matically because, when the observer starts moving, the 
space quantum of the observer moves with it. 

Although the centre of a space quantum can be any¬ 
where in the continuous three-dimensional coordinate 
space, physical space is defined as the space occupied by 
the strings and is therefore quantised. The space quan¬ 
tum of an elementary particle is an extended object and 
Quantum Mechanics also requires that every elementary 
particle is associated with an extended object, namely 
the wavefunction. Perhaps Schrodinger’s 3N-dimensional 
configuration space should not be regarded as an abstract 
mathematical space but should instead be identified with 
the physical space defined by the superposition of N space 
quanta. 
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The inertial mass constraint (IMC) for space quanta 
is introduced in Section [TT] and, in Section III the IMC 
is used to give a physical model for Newtonian gravity. 
When relativistic effects are included and the experimen¬ 
tal constraints are applied, the equations of motion are 
found to be the same as in General Relativity. Inertia 
and kinetic energy are explained in terms of the energies 
in the Machian strings and a Machian formula is given 
for G, the constant of gravitation. Possible solutions to 
the dark energy and dark matter problems are given in 
Sections [TV] and |VJ respectively. 


II. THE INERTIAL MASS CONSTRAINT (IMC) 


Inertial mass is defined as the total energy in the space 
quantum of a massive particle in its instantaneous rest 
frame and the inertial mass constraint (IMC) states that 
the inertial mass of every massive particle is constant. 
The IMC allows an attractive gravitational force to be 
obtained with positive potential energies in the strings, 
as explained in Section III The IMC is also used in Sec¬ 
tion [TV] on ‘dark energy’ and Section [V] on ‘dark matter’. 


III. A MACHIAN MODEL FOR GRAVITY 


The equations of motion for a Machian model of grav¬ 
ity are derived in the Appendix. The model is a direct 
interaction theory, as opposed to a field theory, and is 
defined by the energies in the strings connecting all pairs 
of particles in the observable universe. The energy in the 
strings has the form of positive Newtonian gravitational 
potential energy, with additional terms depending on the 
velocities of the particles. 


A. Inertia and Newton’s laws of motion 


According to Newton’s laws of motion, a free particle 
in an ‘inertial frame’ moves with uniform velocity and a 
force ma is required to give a mass m an acceleration a. 
In General Relativity, acceleration is defined relative to 
the local ‘spacetime’ around the particle and it has been 
suggested that an interaction with the quantum vacuum 
is responsible for the existence of inertial forces mm- 
In the Machian string model, the energy of a particle 
is not concentrated at a point but is instead distributed 
throughout the Machian strings. Inertial frames are iden¬ 
tified with frames moving at uniform velocity relative 
to distant matter, in accordance with Mach’s principle. 
Appendix [B] shows that the kinetic energy of a parti¬ 
cle is part of the gravitational interaction energy in the 
Machian strings. Appendix [D] shows that distant matter 
exerts a force — ma on a particle that has acceleration a 
relative to distant matter, so that an external force ma 
is needed to maintain the acceleration. 


B. Newtonian gravity and the constant of 
gravitation 


In the Machian string model, Newton’s point masses 
and negative gravitational potential energies are both re¬ 
placed by positive gravitational potential energy in the 
strings. In the simplest model of Machian space quanta, 
the direct strings in Figure [I] joining the centres of the 
two masses mi and m 2 , a distance r apart, have en¬ 
ergy +Gm\m 2 /r, where G is the constant of gravitation. 
Due to the action of the IMC, the total energy at the 
centre and in the Machian strings connected to mi is 
mic 2 —Gm\m 2 /v, and similarly for m 2 . The total energy 
in the centres and all the strings of mi and m 2 is there¬ 
fore mic 2 + 7712 c 2 — Gm\in 2 /r, so the interaction energy 
is the same as in Newtonian gravity. The idea that the 
gravitational mass energy is positive, with a compensat¬ 
ing reduction in the energy of each mass, was suggested 
many years ago mm- 

There is a well-known Machian relation for G [5], 
namely 



where Mu and Ru are the mass and radius of the ob¬ 
servable universe, respectively. In the string model, the 
relation ([l]) follows directly from the fact that the iner¬ 
tial mass, me 2 , is the same order of magnitude as the 
total gravitational potential energy in the strings. Since 
the average length of a string is of order Rjj , the total 
gravitational potential energy in the strings connecting a 
mass m to distant matter of total mass Mjj is of order 
GmMu/Ru so me 2 ~ GmMjj/Ru- 


C. Inertial and gravitational mass 


The inertial mass in the space quantum of the i th par¬ 
ticle, mi, is defined as the mass at the centre of the space 
quantum plus the mass in all the strings connected to it. 
The gravitational mass of the i th particle, m*, is defined 
as the mass at the centre plus the mass in the (very large) 
subset of strings connected to distant matter only. When 
all the space quanta are superposed, each pair of space 
quanta has two strings connecting their centres. The to¬ 
tal energy in the two strings is Grhiifij/rij so the energy 
in each of the two strings is Grhirhj/2rij. Since the differ¬ 
ence between inertial mass and gravitational mass is due 
to the mass in the strings connected to nearby matter, 
it follows that inertial and gravitational mass are related 
by the equation 


2 - 2 v—' Gniimn 

m z c = niiC + ) — - 

z -—' Sr,.,- 


( 2 ) 


where S n denotes the set of nearby particles in the system 
under consideration. 
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D. Experimental constraints 

1. The two body and three body problems 

Relativistic corrections for the motion of two gravitat¬ 
ing masses and the acceleration of the Earth and Moon 
towards the Sun in the three-body Earth-Moon-Sun sys¬ 
tem are calculated in Appendix [EJ After applying the 
experimental constraints, the equations of motion are 
found to be the same as in General Relativity. In the 
corresponding Machian string model, approximately 8% 
of the mass of a space quantum is in the Machian strings 
and the remaining 92% of the mass is at the centre. The 
effect of a gravitational field on the precession of a gyro¬ 
scope is considered in Appendix [F] and the usual results 
for the rates of geodetic and gravitomagnetic precession 
are recovered. 


2. Gravitational redshift and the propagation of light 

Gravitational redshift and the propagation of light in a 
gravitational field are discussed in Appendix [G] Photons 
propagate with constant energy in a gravitational field 
because a massive particle propagates with constant en¬ 
ergy and a photon may be thought of as a massive par¬ 
ticle with vanishingly small rest mass. From the relation 
E = hf, it follows that the frequency of a photon prop¬ 
agating in a gravitational field is constant. To explain 
the gravitational redshift, consider an atom of mass m 
at rest at a distance r from a mass M. The energy in 
the strings joining the two masses is GMm/r and the 
IMC implies that the energy at the centre of the atom is 
toc 2 — GMm/r , so the atomic energy levels are reduced 
by a factor of 1 — GM/rc 2 . The frequency of a photon 
corresponding to a given atomic transition is therefore 
lower when the photon is emitted in a gravitational field. 
The frequency remains unchanged as the photon propa¬ 
gates out of the gravitational field, so the photon emerges 
with a lower frequency and therefore appears to be red- 
shifted. Appendix [G] also shows that the speed of light is 
reduced in a gravitational field and Fermat’s principle is 
used to calculate the deflection of a light beam and the 
time delay of radar signals. 


3. Gravitational radiation 

The calculations in Appendix [T] show that the Machian 
strings around a massive body are deformed by the pres¬ 
ence of a nearby mass. When two masses orbit one an¬ 
other, the deformation in the strings varies periodically 
in time and waves in the Machian strings are expected. A 
formula for the rate at which energy is radiated is derived 
in Appendix[j]for the case of circular orbits. Comparison 
with the experimentally verified formula for the rate of 
emission of gravitational radiation in General Relativity 
suggests that gravitational radiation may be attributed 
to the generation of gravitational waves in the strings. 


IV. THE EXPANSION OF THE UNIVERSE 
AND ‘DARK ENERGY’ 

At first sight, the Machian relation @ implies that G 
is a function of time since Mjj and Ru vary with time as 
the Universe expands. It is shown in Appendix [Hj how¬ 
ever, that the relation ([lj is not only consistent with the 
expansion history of the Universe but it can, in fact, be 
used to calculate the expansion history of the Universe. 

On cosmological time scales, the expansion of the Uni¬ 
verse has a significant effect on the number of Machian 
strings in the Universe and on the energy in each string. 
In all space quanta, the number of Machian strings in¬ 
creases with time because the number of particles in the 
observable universe increases as the Universe expands. In 
addition, the mass associated with the potential energy in 
each string has a kinetic energy associated with its reces¬ 
sional velocity. In the early universe, it turns out that the 
kinetic energy of the potential energy is small compared 
with the potential energy so, as the Universe expands, the 
total energy in every Machian string decreases. The IMC 
requires the effects of the increasing number of strings 
and the decreasing energy of each string to balance and 
the result is that the expansion of the Universe deceler¬ 
ates. The expansion rates in the radiation era and the 
matter era are the same as in the conventional ACDM 
model. 



FIG. 2: Scale factor, R(t), in the string model (red) and in 
the ACDM model (green) as a function of the time, t. In both 
models, the present time is at R = 1. The two scale factors 
are almost indistinguishable experimentally, as discussed in 
Appendix [H] 


In the later universe, the kinetic energy in each string 
becomes larger than the potential energy and the IMC 
then requires an accelerating expansion to keep the total 
kinetic energy in the space quantum of a massive parti¬ 
cle constant. The time evolution of the scale factor in 
the string model calculated in Appendix [H] is shown in 
Figure [5] and compared to the time evolution in the con¬ 
ventional ACDM model. 
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V. GALAXY ROTATION CURVES AND ‘DARK 
MATTER’ 


It is known that Newton’s law of gravity, when ap¬ 
plied to the observed mass distribution, does not give 
the correct dynamics on the scale of galaxies and galaxy 
clusters. Either there is a very large amount of still un¬ 
seen dark matter or Newton’s law is incorrect. In mod¬ 
ified Newtonian dynamics (MOND) |T5], the Newtonian 
1 /r 2 force law is replaced by a 1/r force law when the 
Newtonian acceleration is less than a characteristic ac¬ 
celeration, do. It has been shown m that MOND gives 
a very good fit to the observed galaxy rotation curves for 
« 0 » 1.2 x 10 ” 10 m/s 2 . 

In the Machian string model, an acceleration scale of 
order ao is a natural consequence of a contact interaction 
between Machian strings. Consider a test mass in the 
vicinity of a larger mass M. In addition to the Machian 
strings of M, the test mass is immersed within the back¬ 
ground space defined by the Machian strings joining all 
the pairs of particles in the Universe. It is therefore rea¬ 
sonable to suppose that the effect of the Machian strings 
of M on the Machian strings of the test mass depends 
on the ratio of the density of the Machian strings of M 
to the density of background strings. In Appendix [TJ it 
is shown that the ratio of densities is equivalent to the 
ratio of the Newtonian gravitational acceleration to ag. 

Physically, the additional acceleration is due to the 
forces exerted on the centre of the test mass by the 
Machian strings connected to it. Due to the action of 
the IMC, the force between every pair of space quanta is 
attractive so the strings are in tension. In the absence of 
the mass M, the distribution of Machian strings around 
the test mass would be spherically symmetrical and all 
the tension forces on the centre would cancel out. The 
interaction with the Machian strings of M causes the 
distribution of Machian strings around the test mass to 
become asymmetrical, so the forces exerted on the cen¬ 
tres no longer cancel out and the test mass experiences 
an additional acceleration. 


To estimate the magnitude of the additional accelera¬ 
tion, consider the total energy in the Machian strings con¬ 
nected to a test mass m. As noted in Section HID about 
8 % of the rest mass, 0.08mc 2 , is in the Machian strings. 
Dividing the total energy in the Machian strings by the 
average length of a Machian string, Rjj / 2, gives a force 
of about 0.16mc 2 /Rjj. In the limiting case of maximum 
asymmetry, with all the Machian strings pointing in the 
same direction at the centre, an additional acceleration 
of 0.1 6c 2 /Ru would be expected. The acceleration ao is 
approximately 0.4 c 2 /Rjj, where Rjj = 3.6 x 10 26 m is the 
radius of the observable universe found in Appendix [H] 
so the maximum additional acceleration is about 0.4ao- 
In contrast to MOND, which simply postulates a 
change in Newton’s law, the additional acceleration in 
the string model is added to the Newtonian acceleration. 
Thus, although the additional acceleration is significantly 
less than ao, it is still sufficient to explain the MOND phe¬ 
nomenology. A specific model for the interaction between 
Machian strings is considered in Appendix [I] and the max¬ 
imum additional acceleration is found to be about 0 . 2 ao- 


A plot of the Newtonian gravitational acceleration and 
the additional acceleration calculated in Appendix [I] is 
shown in Figure [3j plotted in units of ag. The total ac¬ 



FIG. 3: The Newtonian acceleration (red curve) and the ad¬ 
ditional acceleration due to the interaction between Machian 
strings (green curve) around a point mass M in the Machian 
string model for the model considered in Appendix]!] The ac¬ 
celeration is plotted in units of the MOND acceleration scale, 
ao, and the radial distance r is plotted in units of the radius 
R g , defined as the radius at which the Newtonian gravita¬ 
tional acceleration is equal to ao, so the MOND acceleration 
is equal to 1/r 2 for r < 1 and 1/r for r > 1 . 

celeration in the string model is seen to be very close to 
the MOND acceleration when the acceleration is much 
less than ao, so the string model is expected to give a 
good fit to the observed galaxy rotation curves. 


VI. SUMMARY 

In the Machian string model, direct strings give New¬ 
tonian gravity, Machian strings give inertia and the ex¬ 
pansion of the Universe and dark matter effects are due 
to the interaction between Machian strings. 

A formula for the gravitational potential energy in the 
strings is derived from a relativistic string action in Ap¬ 
pendix |Aj The usual formula for the Newtonian kinetic 
energy of a massive particle can be derived by a suitable 
choice of the fraction, re, of the mass of a space quan¬ 
tum that lies in the strings. It is shown in Appendix [5] 
that the relativistic correction to the kinetic energy and 
all the experimental results for the relativistic two-body 
problem, the Nordtvedt effect and gyroscope precession 
are obtained if the free parameters a±, a 2 , b\ and 62 are 
suitably chosen. The theory gives the standard results for 
gravitational redshift, light bending and radar time delay 
and the electromagnetic interaction between charged par¬ 
ticles can also be included in the formalism, as explained 
in Appendix [Gj A physical model for the observed rate 
of emission of gravitational radiation from binary pulsars 
is discussed in Appendix [TJ 

The expansion history of the universe can be calcu¬ 
lated by considering the total energy in all the Machian 
strings of a space quantum. The usual results for the 
decelerating expansion of the Universe in the radiation 
and matter eras are obtained if the mass fraction re is 
assumed to be constant. By including the effect of the 
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Hubble flow on the energy of a Machian string, a good fit 
to the experimental data is obtained without introducing 
‘dark energy’. 

The interaction between Machian strings gives effects 
usually attributed to ‘dark matter’. The strength of the 
interaction depends on the density of space defined by 
the Machian strings. Since the ratio of the density of 
space around a mass to the background density of space 
is equivalent to the ratio of the Newtonian gravitational 
acceleration to the acceleration scale ao, the MOND phe¬ 
nomenology is explained. 
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Appendix A: Lorentz invariance and the Lagrangian 
for a string 


Let A /i (r, a) be the spacetime coordinates of a classi¬ 
cal string joining the centres of the i th and j th particles, 
where r is a time-like coordinate and a £ [0,1] parame- 
terises the position along the string m- An action for 
the string may be written in the form 

S = f C drda , (Al) 


where the Lagrangian density C is some scalar function of 
the four-vectors dX^/dr and dX^/da. The most general 
form of a Lagrangian density with dimensions r n is 


C = 


\A{8) 

rdX» dX^\ 

to 

{ dr da J 


- B 



2-i n/2 


(A2) 


where A(s) and B(s) are arbitrary functions of s = 
( 9 X m / 9 t) 2 /c 2 . The Nambu-Goto Lagrangian [T5], for 
example, has A = 1, B = s and n = 1. Let the centre 
of mass frame of the string, S cm , be defined as the frame 
in which the velocities of the ends of the string are equal 
and opposite. If the parameter r is chosen to be the time 
coordinate in S cm then 


S — J L cm dT , 

where the Lagrangian in S cm is given by 
Lcm — / L da. 


(A3) 


(A4) 


If x denotes the spatial coordinates in S cm then X M = 
(cr, x) and 

dX^ dX^ 

^ r = (c,x) and ^ = (0,x'). (AS) 


It is desired to find the Lagrangian L in a general 
frame, S , in which the velocities of the i th and j th par¬ 
ticles are v, ; and v ;/ , respectively. In the special case 
B(s) = sA(s), the action is invariant under reparameter- 
isation of the r coordinate of the form r -A t'(t, cr) and r 
can then be chosen to be the time coordinate in the gen¬ 
eral frame S as well as in S cm m- More generally, the 
Lagrangian must first be evaluated in S cm , using (A5), 


and a Lorentz transformation must then be applied to 
find the Lagrangian in S. 

In the Machian string model, the energy in a string of 
length r is proportional to 1/r, so n = —1. Consider a 
straight string in the centre of mass frame. Then x' = £, 
where £ is the instantaneous separation of the two centres 
in S cm , so (dX^/dT)(dX^/da) = -£.x, (dX^/da) 2 = 
—£ 2 and s = 1 — x 2 /c 2 . It is convenient to rewrite the 
Lagrangian density (|A2|) in the form 


£ = /(*) 

M 

s 


f 9(z) 

rdX>‘ 

dX 11 \2 

fdX» \2i 

-1/2 

[ c 2 

K dr 

da ) 

V da ) 



! + <?(*) 


(fx) 2 


- 1/2 


(A6) 

(A7) 


where f(z) and g(z) are arbitrary functions of z = x 2 /c 2 
defined by A(s) = g(z)/f{z) 2 and B(s) = l//(z) 2 and £ 
denotes the unit vector along £. Expanding / and g in 
powers of z, the Lagrangian (A4) has the form 


T f 1 

Lcm = ~z da (1 + a\Z + a2Z 2 + ...) 

4 Jo 


(f xj2 n _i/2 

i + (6 1 + ^ + ...)^r 1 

c z 


(A8) 


for some constants T, ai, a^, &i and b 2 . Let (ct, r) denote 
the coordinates in S. To lowest order, the relative posi¬ 
tions and velocities of the ends of the string are Lorentz 
invariant, so £ may be replaced by ? ij, the unit vector 
along the string in S, and the velocity of the string as a 
function of cr may be written as x = (1/2 — cr)v iJ . Then, 


to (D(v 2 /c 2 ), (A8) gives 


T r a! vfj h (?*. 
Lcm ~ £ i + 12 c 2 24 


+ ... 


(A9) 


If S cm moves with velocity w relative to S then dr = 
dt/'yiw). From (A3), L = L cm /'y(w), so 


L = 


7M£ L 


ai v-j bi (r ij.Vij) 2 


1 + — -ZL _ 

12 c 2 24 


. (A10) 


The quantity l/ 7 (u;)£ must now be rewritten in terms 
of the positions and velocities at a given time in frame 
S. The instantaneous separation £ of the two centres in 
frame S cm can be related to the instantaneous separation 
r ^ in frame S by a modification of the method in m that 
is symmetrical with respect to the two centres. Consider 
the two events A and C in Fig[4j which are simultaneous 
in Scm at time ta and separated by the distance £. In 
S , the events A and C are not simultaneous. A line 
of simultaneity in S is shown in Fig [4] connecting the 
events B and D. The events B and D are defined to be 
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FIG. 4: Diagram showing the worldlines of the i th and j*' 
particles in the centre of mass frame, S C m■ The line AC is 
a line of simultaneity in <S cm and the line BD is a line of 
simultaneity in the frame S. 


For interactions with distant matter, the centre of mass 
frame of rrii and particles of distant matter is effectively 
the rest frame of the distant matter. The calculation of 
relativistic corrections to kinetic energy in Appendix [B] 
requires 0{v 4 /c 4 ) terms proportional to vf/c 4 . The ve¬ 
locity of rrii in S cm is Vjj/(1—Vj.Vj/c 2 ), but the 0(y 2 /c 2 ) 
correction to v^- does not lead to any terms in the La- 
grangian proportional to vf/c 4 and is therefore neglected. 
The velocity of the string as a function of a may there¬ 
fore be written as x = erv^. Similarly, it is sufficient to 


use the approximation £ = in the evaluation of (A81 
and expansion to 0{v 4 /c 4 ) then gives, instead of (A10), 


r 1 1 ,ai ^ _ h '- To -'A/-'’ ,02 v %_ 

7 (w)£ 1 3 c 2 6 c 2 5 c 4 

3 b\ (Yij.Vjj) 4 _ b 2 + «i&i (Vij) 2 {rij-Vjj) 2 ' 

40 c 4 10 c 4 

(A17) 


symmetrical relative to the line of simultaneity in S cm , 
so that t b — ta = ta — te>. The Lorentz transformations 
relating the positions of events B and D in the two frames 
are 


xb = r B + (7 - l)w(w.r B ) - 7W t B (All) 
and X£i = r£) + (7 — l)w(w.r£i) — 7wt£). (A12) 

Since the velocities of the centres are equal and opposite 
in S cm and t b — ta = ~{t b , — tc), it follows that Xg - 
xa = xp — xq. Then 


Evaluating the prefactor (A14) with w = v y instead of 
w = (vj + Vj)/2 gives, instead of (A15), 


7 (w)€ 


, _ l '.i (fij-Vj) 2 

2c 2 2c 2 


with no Q(v 4 /c 4 ) terms proportional to vf/c 4 . 


tuting into (A17) then gives, instead of (A16), 


(A18) 

Substi- 


£ = 


x c = X B - x D = Yij + (7 - l)w(w.ry) 
w(w.r^) 

2c 2 

(r i;j .w ) 2 

' U ‘ r .2 


£ = + 
? 2 = r i 


+ ... 


(A13) 


For interactions with nearby matter, the centre of mass 
velocity is w = (vj + Vj)/2 so 


02 v fi 
5 c 4 


(rij-Vj) 2 01 v ij _ h (rjj.Vjj) 2 
2 c 2 2c 2 3 c 2 6 c 2 

3 b\ (rjj.Vjj) 4 _ b 2 + a 1 b 1 ; v,, S 2 (r, f .v,, ) 2 
40 c 4 10 c 4 


(A19) 


1 

7 (w)£ 



■ 2 , (rij-w) 2 

*3 r .2 


- 1/2 


(A14) 


w 2 

( f b- w ) 2 . 


2 c 2 

2c 2 J 


1 <2 

Vi-Vj (fjj.Vijf 

8c 2 

2c 2 8c 2 




1 


2c 2 



(A15) 


and the Lagrangian (A10) becomes 


- JL h , ( _ lN l ^ _ V *- V J 

" r i:j L V12 8/ c 2 2c 2 

( h J_ In (r,j.v,, j 2 •:i- /f .v,ijf,,.v,j 

V24 8/ c 2 2c 2 


(A16) 


Appendix B: The total Lagrangian and experimental 
constraints 


To define the total Lagrangian for a system of particles 
interacting with each other and with the distant universe, 
it is necessary to define the gravitational mass of the i th 
particle, rhj, for each interaction. For interactions be¬ 
tween a particle in the system and a distant particle, the 
gravitational mass of each particle is equal to the iner¬ 
tial mass, rrii. For interactions between two particles in 
the system, the gravitational masses are defined by equa¬ 
tion ([2]). The total Lagrangian for a system of nearby 
particles interacting with distant matter, using (A9) for 
interactions between pairs of particles within the system 
and (A16) for interactions with distant matter, is there- 






















































If the constant k satisfies the relation 


fore 

*= E 


Gmiifin ' 


V12 8/ 




h , 1 \ (ry.Vi)(%.Vj) 


+ -V 

V 24 8/ 


► E - 

i£S n ,j£S d 


' *3 


V \ 2 

1 - 3 - 

2 c 2 


2 c 2 

fer v i ) 2 

2 c 2 


+ ... 


oi ^ _ 

3 c 2 6 


6 i (Tij.Vjj ) 2 . 02 3bf (jij.Vij) 4 


—f- —— —— 

5 c 4 


40 


b 2 +oi&i (vij) 2 (rij.Vij)2 


(Bl) 


where S n denotes the set of nearby particles in the system 
and Sd denotes the set of particles in the distant universe. 
Let the gravitational potential due to distant matter, k, 
be defined by 


v ' GlJlj 

jes d r b c2 


(B2) 


where r. t j is the distance from a distant particle of mass 
rrij to some nearby mass rrii, and suppose that, in frame 
S , all distant matter particles have velocity v<j. The av¬ 
erage values of (r.y.n ) 2 and (rjj.n ) 4 when summed over 
distant matter are 1/3 and 1/5, respectively, for any fixed 
unit vector n, so (Bl) becomes 


L = y nrriiC 2 
ies n 


9i ; 2 

i _ zzA 

3c 2 


(- 44 + 4 - 


+ E 

'h 
,24 


Grhirhj ' 


/ai l\ V ij v,.v 

1 + (i 2 -sJ^f“ 


V24 8 / 


4 (4- v «j ) 2 


where 


2 c 2 


1 / 3 bl b 2 + aifci \ 

s = d“ 2+ 40- —) 


2 c 2 


(B3) 

(B4) 


From equation ([2]), the gravitational mass rhi is given in 
terms of the inertial mass m,; by 


~ 9 9 V—' Gui^TTlj 

imc = rriiC — > -- 

^ 2 m 


= rriiC 


Gmiirin 


2 r,;.,- 

j£S n ,j^i 3 


- E 


(B5) 


Substitution into (B3) gives 


L = E 


nrriiC 


i£S n 


2v 2 

1 _ fyd 
3c 2 


+ E 

i,jes n ,i<j 


Gmiirij 


1 / &i\ 

sr- 6 J 

i + f-- 1 ) 

V12 8/ 




id 


1 \ V\ 

c 


8 V -^f 

c 4 




2 c 2 


(bi 

+ 4 

(fy.Vy ) 2 

(%-Vi)(f ii .V J ) 

- 

V24 

c 2 

2 c 2 



E 


i,j,keSr,,j^i,k^j 


G 2 mirrijmk 
2 r t jk c 2 


(B 6 ) 


(- - 1) 


(B7) 


and the parameters 02 and 62 in (B4) are such that 

k 5 = — 


(B 8 ) 


then, after dropping the first two terms in (|B 6 j) , which do 
not contribute to the equations of motion, the Lagrangian 
has the form 


L = 


E [ 


ies n 


1 v 4 

^m l v 2 ld + 


E 


Gm, 


1 + /9 


(4 + v 2 ) 


(fij.Vi ) 2 + (r ij .v J -) 2 (r ii .v i )(f ii .v j ) 


+ 


E 


+ V- 
G 2 mimjnik 

rijTjkC 2 


+ • • • , 


(B9) 


i,j,k€S n ,j^i,k^j 

where the parameters /3, 7 , /r, ^ and e are given by 

1 ai 1 b\ 1 

8 ’ 7 = ~~6~ _ 4’ M = ~24 _ 8’ 
61 1 1 


ai 

12 


v = — —- and e = — - 

12 4 2 


(BIO) 


The derivation of the Lagrangian (B91 shows that the 


kinetic energy of a mass is due to the velocity-dependence 
of the energy in the Machian strings connected to it. 

The equations of motion are presented in Appendix [C] 
and relativistic effects are calculated in Appendix [EJ 
The experimental constraints (E 8 ) and (E9) from the 
two-body problem and (E31), (E32) and (E33) from the 


Nordtvedt effect have the unique solution 

3 7 11 

P = 2 ’ T = ~2’ Ai = °’ v = ~2 and 6 = ^2 ' 

(Bll) 


With the parameters in (Bll) and in the special case 


where the reference frame is at rest relative to distant 
matter, the Lagrangian simplifies to 

„ 4 , 


r rl 2 1 1 

L = 2^ l2 miVi + 8 mi ] 

i£S n 


E 

t,j£S n ,z<j 


Grrij 


1 + 


3(4 


2 c 2 

(r (J .v ( !(f,,.v,) 


7vi-vj 

2 c 2 


E 


i,j,keS n ,j^i,k^ij 


2 C 2 

T 2 rriimjmk 
1 kC 2 


(B12) 


which is the same as the Lagrangian derived by Einstein, 
Infeld and Hoffmann m- The Lagrangian ( |A 8 [ ) describ¬ 
ing the motion of the string in its centre of mass frame 
is therefore consistent with experiment if 


39 


ai = 


and bi = —3 . 


(B13) 
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Equation (B7) then gives 


k = 


40 


and (B 8 ) requires a 2 and 62 to satisfy 


a 2 


6 


251 

120 


(B14) 


(B15) 


Equation (B14) implies that the total energy in the 


Machian strings of a space quantum is 

E 


Gmiirin 3 2 

- = m miC > 

Tin 40 

jes d 13 


(B16) 


so that approximately 8 % of the energy of a space quan¬ 
tum is in the Machian strings and the remaining 92% is 
at the centre. 


Appendix C: Equations of motion 


gives the equation of motion for m,, 


Variation of the Lagrangian (|B9|) with respect to r.; 

(Cl) 


dL d / dL\ 
dr i dt V dvj) 


For the mass mi, the equation of motion is 


ies„ 


r iV 


E 


3 ^ (r H .Vi ) 2 + (rij.v ,;) 2 + 3 ^ (rH.vi)(r H .Vi) 
c 2 c 2 

Gmirrii r (rij.v Qvi + (r^.v^v* 


rue 


ru 


+ v 


(fii.Vj)vi + (r li .v 1 )v i 


+ 


chq 


E 


r u 

G 2 mim-jmk 


-' TinTnkC 


-dt) mi 


(-S) 


E 

ieSn 


2 c 2 
Gmimj 


Vld 


ruc A 


2/3v! + yvi + 2/i(r li .v 1 )r li 

+ | = 0. (C2) 

Appendix D: Inertia and Newton’s laws of motion 

Suppose that the masses m.\ and m 2 are in a reference 
frame accelerating with respect to distant matter with 
acceleration a. In the accelerated frame, distant matter 


has an acceleration —a, so f d = —a. Equation (C2 1 then 
has the form 


Gmim.2 

rriiYi = —mi& -—r 


(Dl) 


Equation (Dl) defines the total force acting on the mass 


mi in an accelerated reference frame. When m 2 = 0, 
the additional force that has to be added to the right 
hand side to give iq = 0 is mia. A force F = mia is 
therefore needed to maintain an acceleration a relative to 
distant matter, which is Newton’s second law of motion. 
Newton’s first law of motion, that a particle moves with 
fi = 0 in an inertial frame when there are no external 
forces, corresponds to the special case a = 0. 


Appendix E: Relativistic gravity 

1. The two body problem 

Consider the motion of two masses, mi and m 2 , in an 
inertial fra me. Setting Td = 0 and = 0, the equation 
of motion (C2) for mi is, to first order in c 2 /c 2 , 


ri 


1 (■■ v( . 0 vi.iq 

r 


2 \ c 
G«2 


i +/J q+^ +7 ^ 

(f.vi ) 2 + (r.v 2 ) 2 D (r.vi)(f.v 2 ) 

3/i- 2 -b - 2 - 


Gm 2 r (r.vi)vi + (r.v 2 )v 2 


rc* 


Gmo 


(r.v 2 )vi + (r.vi)v 2 


+ r 3 c 2 (r.v) 2/3vi + qv 2 + 6/t(r.Vi)r + 3^(r.v 2 )r 

Gm 2 b... .. „ (v.vi)_ „ (r.vi) 

-— 2 /frq 4- yr 2 + 2 /u --r + 2 / 1 ---v 


rc* 


_ _ (r.V 2 ) (v.v,V _ _ 

+ 2 /r(r.ri)r + v- - -v + v- --r + r^(r.r 2 )r 


2eG 2 


( mirri 2 + m 2 )r . 


(El) 


After substituting the Newtonian appro xima tions iq = 
-Gm 2 r/r 3 and r 2 = Gmir/r 3 , equation (El) reduces to 


iq — A 1 r + B [Vj + B 2 V 2 


(E2) 


with 

Gm 2 G 2 m 2 ' 

Ar — ———h 

Gm 2 r 


(2/3 + 2 /i)m 2 - (7 + v)mi 


- /3(v 2 + vl) - 7 V 1 .V 2 

- 3/i {(r.vi) 2 + (r.v 2 ) 2 } + 6/t(r.v)(r.Vi) - 2/iv.vi 

- 3ii(r.V!)(r.v 2 ) - v(v.v 2 ) + 3zi(r.v)(r.v 2 ) 

G 1 TI 2 2 2eG 


+ ^2 U 1 - 73^( TO l TO2 + W 2) > 


Bi = 


and B 2 = 


Gm 2 r 


Gm 2 ' 


r.vi + 2/Jr.v 
(v + 2/i)(r.Vi + r.v 2 ) + yr.v 


(E3) 


The equation of motion for r in the centre of mass 
frame may be found by subtracting the corresponding 
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equation for m2 and substituting Vi = TO2V/M and 
V2 = —miv/Af, where M = m 1 + m2, which gives 


with A = — 


r = Ay + B (r.v)v 
GM G 2 


(E4) 

(2 (3 2/r — 2e)(m 1 -t- m2) 

— 2(7 + z/ + 2e)mim-2 


G 

Mr 2 c 2 - 


{3/z(m 2 + m 2 ) — 3z/mim 2 }(r.v) 2 
|(/3 + 2/x - (mf + to 2 ) 


- (21/ + 7- ^ )mim 2 }> f 2 


and B = 


G 

Mr 2 c 2 - 


(1 + 2/3) (to 2 + m|) 

— (1 + 27)mim 2 . (E5) 


Introducing pol ar coordinates for r, the angular com¬ 
ponent of (E4) may be integrated to give r 2 6 = 
hexp[f' B dr' J, for some constant h, and the radial com¬ 
ponent becomes 


d 2 u 

d(P 


+ u = - 


h 2 u 2 

GM G 2 i 

+ 


o-2 S r Bdr' 


h 2 


(2 -(- 2,(3 — 2/i -(- 2e)(mJ -t- m^) 


h 2 c 2 . 

— (2 + 27 — 2v — 4e)mim2 
3G[^mim2 — /z(m 2 


G 


Mh 2 c 2 L 


Mh 2 c 2 

1' 


m l)\ ■2 


(/3 + 2/r- -)(to? + to|) 


- (7 + 2i/ - -jmim 2 


(E6) 


where u = 1/r. After using the energy conservation equa¬ 
tion u 2 = 2GMu + 2EM/mirri2 to substitute for v 2 , the 
rate of periastron precession is found from the coefficient 
of u on the right hand side to be 


A d 


7tG 2 

h 2 c 2 


(1 -(- 4/3 -f- 2p -(- 2e)(m 2 -f- to*]) 


— (1 + 4y + 2v — 4e)miTO2 


(E7) 


per orbit. The prediction of General Relativity is Ad = 
67 tG 2 M 2 /h 2 c 2 , so agreement with the experimental data 
for the precession of the perihelion of Mercury and the 
periastron of a binary pulsar requires 

4 (3 + 2p + 2e = 5 (E8) 

and 4y + 2v — 4e = —13. (E9) 


2. The Nordtvedt effect 

The absence of a Nordtvedt effect in the Earth-Moon- 
Sun system 20] implies that the accelerations of the 
Earth and Moon towards the Sun are the same, so that 


the acceleration of a composite body towards the Sun is 
independent of its gravitational self-interaction energy. 
For a composite body consisting of two point masses, m* 
and rrij . separated by a distance of r t j , the acceleration 
of the centre of mass towards a massive external body 
must therefore be independent of Gmiirij/rij. 

Consider the motion of the two masses nii and rrij in 
the presence of an external mass m e . The acceleration of 
the centre of mass is given by 


a = 


mAi + rrijVj 


M 


rrii .. 

M ri 




(E10) 


where M = rrij + rrij . The acceleration f, of m; in the 
prese nce of rrij and m e may be found by extending equa¬ 
tion (El) to the case of two nearby masses. The required 


generalisation of (Ell may be written in the form 




Grrij 
TijC 2 l 

Gm, 


2 eG 2 


2/3?* + 7 fj + 2n(r ij .r i )r ij + v(rijXj)r i: 
2/3r i + 7? e + 2/x(r ie .r i )r ie + v(r ie .r e )r 


uiirrij + m 2 rrijm e / 1 1 


X, 


13 

m l m e + ml 

r^ 
ie 

+ Xi e , 


\r ie r je y 


/ I 1 m 1 'i 

— + — r *e 

\Vin r^pJx JJ 

(Ell) 


' je' 


where X, 7 denotes the remaining acceleration- 
independent terms of the right hand side of (Ell 
corresponding to the interaction of rrii and rrij , i.e. 


X« = - 


Grrij 


1 + 


{Vi+V 2 ) Vj.Vj 

f -2 — -1" 7- 

C z c z 


+ 3/i 


(Tij.Vi ) 2 + (r jj.Vj ) 2 


Grrij 


(f+ (r 


2^ 13 




Grrij . . r _ __ 

+ -3 —2/3v* +7Vj + 6/z(r ii .v i )r ij - 

r b' c L 


Grrij r (v.^-.Vj), 

-1 2 

T ijC T ij 


2/i 


+ 3is(rjj . Vj )r, 

(rjj.Vj) 


T ^ -l 


■ (E12) 


The acceleration-dependent terms in (Ell) may be eval- 
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uated by substituting the Newtonian approximations 


Gm,- 


r i = -- 


Gnrii „ 


2 


and 


Grrii 


Gm e 


Gm e ^ 
r 2 r l e > 


r R = 


9 Me 


r,: e + 


Grrii 


(E13) 




Substituting the Newtonian approximation for r* into the 
first line of (Ell) gives 


f< = Y 


Y, 


2eG 2 crminrij + m? mjm e / 1 1 

c 2 U rfj + r« a V 


' ie ' je 


miTYle + m e 


+ S^(± + J_)| r 1}, (E14 ) 

Vi* 3 \r 7 „- Va p J\ \) 


where 


Y - x | Gm i 
lJ ~ + 2r 2 -c 2 


{h 2 bj + (vi.rjj)vj} 


Grtij r ... .. ,, „ .. 

-y 2 prj + 717 + 2 //(r,,.r, ir,, + 

TijC L J 

(E15) 

To evaluate Y ij, it may be noted that if the parts of 
r, : a nd r, propo rtion al to ry in (E13) are substituted 
into (E15) then (E15) has the same form 


A i*i 


B{ Vi 


B 3 2 Vj 


(E16) 


as given in (E2), where A \, B{ and B J 2 are of the form 
corresponding to (E3) for the interaction of rrii and rrij , 

i.e. 


A i — ~ 2 

r l 


Grrii G 2 m.i ' 


+ 


r 3 „2 


(2/3 + 2 n)mj — (7 + v)rrii 


+ 


Gnrii 


~~ P( v i + Vj) — 7 v i- 


- 3/i {(r,j.v ,) 2 + (r^.Vj) 2 } + 6 /i(r ij .v ii )(r ij .v i ) 

- 2 fJ.Vij.Vi - 3^(fjj.Vj)(fjj.Vj) - ^(vij.Vj) 


+ 3z/(r y .Vij ) (r jj .Vj) 


B l = S5 ? b-Vi + 2^?ij.Vy 


Gnrij 2 


2^c 2 




. Grrij r 

and ^2 = -2-y, 


+ 2 /x)(?jj.Vj +r ii .v j ) 
+ Tf«-vy ■ 


(E17) 


Adding the additional contributions to Yy from the 
parts of Vi and Vj in (E13) proportional to r, e and r^e, 
respectively, gives 


Y./ — Ty 


13( Vi + 


+ 


Gm.j f Gm e 




V 


2/3r ie + 2n(r l j.r ie )r z 


+ 


Gm e ' 

~7 r ~ 

je 


7 ? ie +v{v i jXje)v i j }. (E18) 


Similarly, 

Y ie = Afr ie + Bfvi + B 2 Vj 
Grrip [ Gm 


TieC “ 


2/3?y + 2/r(r le .Ty)r le 


Gmj ' 


(E19) 


To lowest order, the centre of mass acceleration (E10) 


is equal to — Gm e r/r 2 where here, and in the remain¬ 
der of this section, r denotes the position vector from 
m e to the centre of mass of rrii and rrij. The Nordtvedt 
effect concerns corrections to Gme/r 2 that depend on 
the binding energy, — Gmirrij/rij , of rrij a nd rrij. Con¬ 
sider, therefore, the contributions to (E 10 ) proportional 
to Gm e /r 2 . The contribution to (E10) from the 0(G 2 ) 
terms in (E14) that are proportional to m e is 


2eG 2 rrii r 

rrijm e 

f 1 

^ 1 M 

Me 2 l 

[ r ij3 

''he 

H 

^ 1 
CD | 

'm l m e 

rrijm e 

1 

+ 1 ^1 

- 

J- 

Ge 3 

Kr ij 

r je )\ 

+1 





■]}• 


(E20) 


and the ry — dependent part proportional to Gm e /r 2 is 
2eG 2 " 


miinjm e 


1 ' *7° 

4eG 2 


( r -f + r -f) 

V r - r J 

le je 


mirrijm e 


Mrijr 2 c 2 


r. 


(E21) 


The contribution to Y .7 proportional to Gm e /r 2 is, 


from (E18), 

Gtoo-to 


j 11 be 


(2/3 + y)r + (2/i + !/)(%. r)r„ 


and the corresponding contribution to (E10) is 
2 Gmimjm e ' 


MrijV 2 c 2 


(2/3 + 7 )? + (2/j, + v)(rij .?)?» 


(E22) 


(E23) 


The expression (E19) for Y, e contains contributions pro¬ 
portional to Gm e /r 2 from A\ r i e +Bf Vi+B 2 Vj and, from 
the other terms, contributions proportional to Gm e /r. 
The contributions proportional to Gm e /r 2 are 


Gm e Gm e 

o "f - n o 


-Hvf+Ve) 

- 7 Vj.Ve - 3 (1 {(r.Vi ) 2 + (r.v e ) 2 } + 6 /x(r.v ie )(r.Vj) 

- 2 /JLV ie .Vi - 3^(r.Vj)(f.V e ) - v(v ie .v e ) 

Gm e 2 ' 
2 r 2 c 2 Vl >T 


Gm e 


Grrip r 


+ 3z/(r.v ie )(r.v e ) 

+ 2/ir.v ; 

(v + 2 /i)(r.v.j + r.v e ) -(- 7 ?.^ 


(E24) 
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Suppose, for simplicity, that the masses mt and rrij are 
in a circular orbit. Then v 2 is constant and equal to 
Gm^/rij. The time averaged valuesoTr.Vi and v.i.v e are 


zero, so the time-averaged part of (E24) that depends on 


Gm e ^ G 2 m 2 1 m e M 

~^~ r + Mr ijr 2 c 2 [2 ~~ M 

+ 3/i((r .Vi) 2 ) 

G 2 m 2 m e 

3 r(l + 2/3)((r.Vj)vi) 

up-- 


(E25) 


where the angled brackets denote the time average. After 
replacing Vj by Vy, the contribution of (E25) to (ElOl is 


Gnip 


G 2 m,'miTO e r 1 

o ~ P - 


3^((r.v ij ) 2 ) 


rtiirrijme . ^ . 

-|^-(l + 2/3)((r.Vy)vy), 


Mrijr^c 


(E26) 


so the contribution toY.j e proportional to Gm e /r 2 gives 
a contribution to (E10) proportional to Gm e /r 2 . The 
contribution to Y.; e proportional to Gm e /r is 


G 2 rrijm e ", 

r ij r ieC 2 


2 Aty T 2/i(rj e .r. I? )r,; ( 


which gives a contribution to (ElOl of 
G 2 m e mimj 


Mr 2 jC 2 


2/3 (— — ) 

L W ie 


2h{ 


(r le .r.y) 


(rje.r^-). 


■}] 


G2m - m>m ' [/(r.) - /(r,) 


Mr 2 -c 2 

L J 


where /(x) = 2/3 7 


+ 2p 


(x - X e ) .Yjj 


(x-x e ). 


(E27) 


(E28) 


With the approximation /(r^) — /(r 7 ) = (r^ — r 7 - ).V/ (r,;), 
the time-averaged contribution to (ElOl from (E27) be¬ 
comes 


m e rriimj 


2/tr + (2/i - 2/3){(r.r i j)r ij ) 
2\ 


6/t((r.r ii ) 


(E29) 


so the contribution to Y 7e pr oportional to Gm e /r also 
gives a contribution to (ElOl proportional to Gm e /r 2 . 


The time averages of (r.ry)“ and (r.v,j) 2 are equal 
to (sin 2 0)/2, where 9 is the angle between r and the 
normal to the orbital plane of rrii and rrij, and the 
time averages of (r.r, :j - )r,j and (r.Vy)Vy are equal to 
(sin 2 Or + sin0cos0n)/2, where n is a unit vector per¬ 
pendicular to r in the plane containing ? and the normal 


to the orbital plane. The sum of (E21), (E23), (E26l 
and (E29) gives, for the total dependent centre of 


mass acceleration proportional to Gm e /r 2 
Gm u 

+ Mrijr 2 c 2 IA2 +d/: 

3 


— 4e 


[{- +3^ + 2 7 - 

— ^/isin 2 6>j r 


1 


+ -(1 + 6/r + 2zz)(sin 2 9 r + sin 0 cos 9 n) 


(E30) 


The constraints on the parameters needed to ensure the 
absence of the Nordtvedt effect are therefore 


3/3 + 2y - 4e = -i , 

H = 0 


and 


2 ' 


(E31) 

(E32) 

(E33) 


Appendix F: The precession of gyroscopes 

The precession of a gyroscope in orbit around 
the Earth may be calculated using the EIH La- 
grangian ( |B12 ). The method is described in (21] and 
explicit calculations are given here for completeness. 


1. Geodetic precession 

Consider a gyroscope rotating with angular velocity u> 
and moving with velocity V relative to the Earth. The 
Lagrangian for the gyr oscop e may be obtained by inte¬ 
grating the Lagrangian (G1) for a point particle over the 


mass distribution of the gyroscope. A point with position 
vector x relative to the centre of mass of the gyroscope 
has velocity v(x.) = V + ui x x relative to the Earth. 
The geodetic precession comes from the change in the 
Lagrangian that is proportional to the angular velocity 
of the gyroscope and the mass M of the Earth. After 
substituting into (Gil, the last term in (Gl) becomes 


3 GM 
2 c 2 


a P( X M X ) 2 




(FI) 


where p(x) is the mass density of the gyroscope and r is 
the position vector of the centre of the gyroscope relative 
to the centre of the Earth. The required change in the 
Lagrangian is therefore 


6L = 


3 GM f ^ p(x)V.(w x x) 
c 2 1 X 1 ' 1 

3 GM 


VitijkLOj J d 3 xp(x)x k ^ ^ 


(F2) 
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The first term vanishes by definition of the centre of mass 
and, if the gyroscope is assumed to be a sphere, 

J d 3 x p{~x)xiXj = , (F3) 

where / = ^ J d 3 xp(x)x 2 (F4) 


is the moment of inertia. Equation (F2) then gives 


_ 3 IGM 

^ o o * j 


2 r 3 c 2 

3 GM 
2 r 3 c 2 


S.(r x V), 


(F5) 


where S is the angular momentum of the gyroscope. The 
additional term (F5) in the Lagrangian gives a precession 


of the gyroscope with angular velocity 
n __ 3 GM 

^ 9eo ~ 2 r 3 c 2 ( r X V )’ 
in agreement with the experimental value. 


(F6) 


2. Gravitomagnetic precession 

When the rotation of the Earth is taken into account 
there is an additional contribution to the rate of gyro¬ 
scope precession, known as gravitomagnetic precession, 
associated with the angular momentum of the Earth. 
Let the angular velocity of the Earth be u and a con¬ 
sider mass element m\ in the gyroscope with position 
vector x relative to the centre of the gyroscope and a 
mass element m 2 in the Earth with position vector xe 
relative to the centre of the Earth. The velocities of 
the two mass elements are i>i(x) = V + u x x and 
^2 (x_e) = we x xe, respectively, and their relative po¬ 
sition vector is r 12 = r + x- xg. The terms in the La¬ 


grangian (B12) that are proportional to and 

quadratic in the velocities Vi and v 2 are 


GtoiTO 2 

ri 2 c 2 


3 7 1 

~(vf + vl) - -Vi.v 2 - -(ri 2 .vi)(f 12 .v 2 ) 


(F7) 


After substituting for Vi and v 2 and integrating over the 
mass densities p(x) and pe{xe), the terms containing 
factors of both w and u> E are 

8L = -^2 J J d 3 xd a x E p{x)p E (x E ) 

, (u> x x).(u E x x E ) 

ri2 

r 12 .(u> x x)r 12 .(w £ x x E ) 


12 


= f d 3 x p(x)I(x) 


(F8) 


where, after defining R = r + x, the integral 7(x) may 
be written in the form 


J(x) = J d 3 x E Pe(*e) 


,{oj x x) . (uj E x x E ) 


r 12 


R.(w x x)R.(ug x xe) xg.(w x x)R.(wg x x E ) 

(F9) 


12 


' 12 


The integral (F91 may be evaluated as a power series in 


1 /R using the expansions 


1 R.xe 3(R.xe) 2 — ■ 


n 2 R 


R 2 


2 R 3 


, 1 1 3R.xe 

and 1 


12 


R 3 


R 4 


+ ... 

(F10) 


Assuming the Earth to be spherical, the integral of a term 
containing two factors of xg may be evaluated using the 


formula analagous to (F3) and the integral of a term 
containing an odd number of factors of xg is zero. The 
expansion of (|F9|) reduces to 


/( x ) = Me 


R.[(w x x) x we] 

r: 2 


(Fll) 


to order 1/i? 2 , and substituting back into (F8) then gives 


sl = (F12) 

c z I R z 


The integral ( |F12 ) may similarly be evaluated in powers 
of 1 /r using the expansion 


R 

R? r 


r x — 3r(r.x) 

~2 ' W3 


(F13) 


which gives 

SL= 2 -^ 


J d 3 xp(x) |[w(o7e.x) -x(o7.a > E )] 

,[x-3r(r.x)]}. (F14) 


After applying (F3), equation (F14) reduces to 


SL = ^--^-[bj.uiE — 3(f.w)(r.wg)] 


G 


S.[J — 3r(f.J)], 


(F15) 


where S is the angular momentum of the gyroscope and 
J is the angular momentum of the Earth. The rate of 
gravitomagnetic precession is therefore 


Llnm. — 


G 


“am - r 3 c 2 


[3r(r.J) - J], 


(F16) 


in agreement with the standard result. 
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Appendix G: The propagation of light 


The variational principle therefore becomes 


Consider a test mass m moving in the gravitational 
field of a mass M. After setting t he vel ocity of the mass 
M equal to zero, the Lagrangian (B12| for the mass m, 
in a frame at rest relative to distant matter is 


L = —mv H— 


1 mv i GmM / 3v 2 \ 

- -(l + ^), (Gl) 


8 e 2 


to lowest order in GM/rc 2 . Aft er subtracting the con¬ 
stant me 2 , the Lagrangian (Gl) may be written in the 
form 


L = — 


1 - 


c(r) 2 . 


1/2 


r)c(r) 2 


(G2) 


for suitable functions m{r) and c(r), which is a gener¬ 
alisation of the Lagrangian —me 2 1 7 for a free particle. 
Expanding the gamma factor in ( |G2| ) and comparing 
with (Gl) gives 


/ \ / \ 2 2 L GM \ 

myr)c{r) = me I 1-— I 

and m(r) = m^l + 


rc 2 7 
3 GM^ 


re 


(G3) 

(G4) 


s[^=0, (G10) 

J c(r) 

where the photon frequency is constant along the path, 
which is Fermat’s principle of least time. 

A photon may also be considered as a localised 
quantum of the electromagnetic field since, as is well 
known [23] . Fermat’s principle is the geometric optics 
limit of Maxwell’s equations for an electromagnetic wave 
in a dielectric medium in which the speed of light is 
c(r). However, the propagation of electromagnetic waves 
through the entire coordinate space is not consistent with 
the string model, since physical space is restricted to the 
strings. For consistency, electrodynamics should be for¬ 
mulated as a direct action theory, in which charged par¬ 
ticles interact directly with each other via the strings and 
there are no independent degrees of freedom associated 
with the electromagnetic field. Classical electrodynamics 
can in fact be reformulated as a direct action theory by 
eliminating the electromagnetic field variables from the 
action [21]. For two particles, with charges q\ and <72, the 
resulting Lagrangian to 0(v 2 /c 2 ) is 


where (|G3|) comes from equating the terms independent 


of v and ( |G i| ) comes fro m eq uati ng th e terms propor¬ 
tional to v 2 . Combining (G3) and(G4| gives, to lowest 
order in GM/rc 2 , 


z{r) = c( 


1 2GM ^ 


rc* 


(G5) 


The momentum corresponding to the Lagrangian 


is 


p = m(r)v 


and the Hamiltonian is 


1 - 


c(r) 2 


-1/2 


H = 


1 - 


c(r) 2 . 


-1/2 


m(r)c(r) 2 


(G6) 


(G7) 


A photon may be considered as the zero-mass limit 
of a massive particle. If the limit is taken so that the 


energy remains constant then equation (G7) implies that 
the velocity v tends to c(r), so c(r) is identified as the 
speed of light in a gravitational field. The variational 
principle needed to calculate the path of a photon may be 
derived by considering the zero-mass limit of Hamilton’s 
principle. Since the speed of the photon is given by c(r) 
on the varied path as well as on the actual path, the times 
at the endpoints are not fixed during the variation. It is 
therefore necessary to rewrite Hamilton’s principle as a 
principle of least action [22] , which states that 


L = 


£[ 


1 


1 


,, 4 -, 


-m lVi + -mj-j 


<M2 

r 


glg2 
2rc 2 - 


vi-v 2 + (r.V!)(r.v 2 ) , (Gil) 


which is known as the Darwin Lagrangian. The La- 
grangia n (|G11 ) can be included in the formalism of Ap¬ 
pendix ALby taking T = —qiqj/&'Keo, a\ = 3/2 and 
b\ = —3 in the Lagrangian (A16) . 


1. Gravitational redshift 


Equation (G3) shows that the mass-energy of a test 


mass at a distance r from a mass M is multiplied by 
a factor 1 — GM/rc 2 . If the test mass is an atom, 
all the energy levels of the atom are multiplied by the 
same factor and the frequency of light emitted by the 
atom for a given atomic transition therefore has the form 
/ = /o(l — GM/rc 2 ), where fo is the frequency of light 
for to the same atomic transition in the absence of the 
mass M. A photon emitted by an atom at a distance r 
from the mass M therefore propagates with a constant 
frequency /o(l — GM/rc 2 ). Compared to the frequency 
corresponding to the same transition in an atom receiv¬ 
ing the photon at a distance r + Ar, the frequency of the 
photon is smaller by an amount A/, given by 


6 


p.dx = 0 


(G8) 


when the path is varied so as to keep the endpoints fixed 
and the energy, H, constant. Now p.dx = p.vds and, in 
the limit v —> c(r), equation (G6) gives 


Hv H 

c(r) 2 c(r) ' 


(G9) 


A/ _ GM GM 
fo rc 2 (r + A r)c 2 

GM Ar 

^ 9 5 

rc z r 


(G12) 


in agreement with the usual formula. The explanation 
for the gravitational redshift is essentially the same as in 
the dielectric model for gravity considered by Dicke [25]. 
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2. The bending of light 


If r) denotes an arbitrary parameterisation of the pho¬ 
ton path the n, in p lane polar coordinates, the path length 
in equation (G10) is given by ds = Vr 2 + r 2 9 2 drj. The 


corresponding Lagrangian, L, is independent of 0 and 
dL/d6 = constant gives 


r 2 6 


= A(l — 


2GM \ 


_|_ r 2 Q 2 V re 2 )' 
where A is a constant. Thus, to lowest order in GM/rc 2 , 


to lowest order in GM/rc 2 , where s is the path length 
measured from the point of closest approach to the Sun. 
If 6 is the angle between the direction of propagation of 
the radar beam and the direction of the Sun then s = 
d cot 9 and r = d cosec 9 and the first integral becomes 


4 GM r r ° ds 
c 3 Jo r 



In (cot 6 + cosec 0) 


cot 1 (r e /d ) 
7t/2 


4 GM 



(G17) 


r dr 
\d9 
/du 
\d9 


2 

2 



A 2 

1 

A 2 


('+ 

( 1+ 


4GM\ 
rc 2 / 
4 GM\ 


(G13) 


where u = 1/r. The zeroth order solution is u = Uq = 
A _1 sin(0 — 9q), so the constant A is identified as the 
distance of closest approach. The first order solution is 
u — uq T u\ where u\ satisfies 


The total time delay of the radar beam is therefore 


St 


4 'GM ^ / 4r e r s \ 
c 3 V d 2 ) ’ 


in agreement with (G15|. 


(G18) 


Appendix H: The expansion history of the Universe 


^ cos(6» - e 0 ) + til sin(0 - 9 0 ) = sin(0 - 0 O ). 

Dividing by cos 2 (9 — 9 o) and integrating gives 
1 2 GM 

u = — sin(0 - 9 0 ) + ^ 2 2 [1 + acos(0 - 0 O )], (G14) 

where a is a constant of integration. The values of 9 — 9q 
corresponding to u = 0 are, approximately, —2 GM/Ac 2 
and 7T + 2 GM/Ac 2 . The angle of deflection of the light 
ray is A 9 = 4 GM/Ac 2 , which is twice the Newtonian 
value, in agreement with observation. 


3. Radar time delay 


Experiments have shown that the time taken for a 
radar beam to travel from the Earth to a satellite is 
slightly increased if the beam passes close to the Sun. 
Let the distance of closest approach of the radar beam to 
the Sun be d and let the distances of the Earth and the 
satellite from the point of closest approach be r e and r s , 
respectively. If the beam travels to the satellite and is 
reflected back to Earth, the time delay due to the pres¬ 
ence of the Sun is found, experimentally, to agree with 
the formula 


St: 


AGM 


ln ( 4 ^ i ) ( d< ^ r e,r s ), 


(G15) 


where M is the mass of the Sun. The total time delay 
due to the red uction in the speed of light from c to c(r) 
is given, from (G10), by 


St = 


1 1 

-c(r) c 

r r ° 2 GM 


ds 


ds+- 


2 GM 


ds, (G16) 


r o 


The expansion of the universe may be calculated by 
taking into account the dependence of the string energies 
on the Hubble parameter. It is reasonable to suppose 
that, if the peculiar velocities of the particles are ignored, 
the energy in the string joining a particle of mass mt to 
a distant particle of mass rrij has the form 


Grriirrij 



(HI) 


for some function /, where H is the Hubble parameter 
and r is the length of the string. For the tests of the 
theory of gravity considered in Appendix [Ej the time 
scales are much smaller than the time scale for the ex¬ 
pansion of the Universe and the effect of the function / 
is simply to rescale the masses rrij of distant matter by 
a constant factor. The calculations in the previous Ap¬ 
pendices therefore remain valid provided the factor / is 
understood to be in cluded in the masses rrij. In particu¬ 
lar, equation (B16) for the total energy in the strings of 


a space quantum is replaced by 


E Grriimj 

T' ■ ■ 


j&Sd 


K 


Hr 


ij \ 2 

1 = nrrijC . 


(H2) 


where k = 3/40. 

To calculate the function /, consider the kinetic energy 
in the strings associated with the Hubble flow. Since the 
Hubble velocity may exceed the speed of light, the kinetic 
energy will be calculated using the Newtonian formula. 
The recessional velocity of a point on the string at a 
distance x from rrij is Hx and the energy of an element 
of string of length dx is Grriimjdx/r 2 j, assuming the mass 
of a string is distributed uniformly along its length, so the 
total kinetic energy in the string is 



1 Gminij 


r 2 . „2 


{Hx) 2 dx 


Gmjin.j tHrjj \ 2 ^ 

6r •ij V c ' 


































16 


and equation (H2) has the form 

Gmi.m , r \ (Hri 


E 

j^Sd 


1+ 6 V c 


= KTYliC 


(H4) 


where y = R/Rq and Rq is the value of the scale factor 
at the present time. It is convenient to define the di¬ 
mensionless time, r, by r = Hot and the dimensionless 
conformal time, 77 , by the equation 


The precise formula for the kinetic energy in the strings 
connected to distant matter may be more complicated 
since the nature of the instantaneous interaction trans¬ 
mitted by the strings is not yet understood. It is therefore 
of interest to consider a more general function / of the 
form f(x) = 1 + ax + bx 2 . Then 


V(t) = H 0 f . (H12) 

Jo y(s) 


Then HqRu/c = yr/ and HR/r/c = 7777 , where the dot 
deno tes d ifferentiation with respect to r. Substituting 
into (Hll) and completing the square gives 


E 

jes d 


Gmirrij 



= KUliC 


2 


(H5) 



4k 

W b ’ 
(H13) 


The sums in (H5) may be evaluated in the continuum 
approximation, which gives 


jes d 


rRu 


P_ dV= IMv 


r 

rRu 


2 Rr 


and 


jes d 


iji l0 


rnu o 

/ prdV = -MuRu , (H 6 ) 
Jo 4 


where M/j is the mass of the observable Universe, R/j is 
the radius of the observable Universe and p is the average 
mass density, and 1H5|) then reduces to the equation 


where y eq = Q r /£l b is the value of y at which the baryon 
and radiation densities are equal, which may be solved 
for y to give 


~ 2k y / VeqY 1 « 2 _ JM 1/2 _ _a 
-3bQ b r] 2 \ y ) 36 b 2 2b. 6 b 

(H14) 


To integrate numerically over several orders of magni¬ 
tude in time, it is convenient to change t he ind ependent 
variable from r to u = In 77 . Equation (H14) then be¬ 
comes 


3 GM V 
2 Ru 


1 + y 


2a ^ HR/j ^ 


2&( 


HRu 


(H7) 


Equation (H7) may be used to calculate the expansion 


history of the universe if it is assumed that the constant 
k, which represents the fraction of the energy of a space 
quantum in the strings, is time-independent. The expan¬ 
sion history is given by the scale factor, R(t), is related 
to the radius R/j by the usual formula 


Ru(t) = R(t) f 
Jo 


R(s) ■ 


(H 8 ) 


The total mass M/j in the observable universe is given 

by 


4:7T q 

M/j = — R/j{Pb + pr) , 


(H9) 


1 dy 

y du 


2k y r Vegy 1 a 2 
-3bttb e 2u \ y ) 36 b 2 


1 n 1/2 
2 b. 


6 b 

(H15) 


The value of fl r is taken to be 8.2 x 10 -5 and is treated 
as a free parameter. The integration begins at the current 
epoch, y = 1, and continues back to the singularity at 
y = 0. The value of 77 at the current epoch, 77 = 770 , is 
chosen so that the Hubble parameter is equal to Hp at 
the present time, i.e. so that y/y = 1. From (H13), it 
may be seen that the required value of 770 satisfies the 
equation 


Vo ( 1 + 2/eg)( 


1 + — 


2 a 
3“ 



4k 

m~ b 7 

(Hi 6 ) 


where p b and p r are the densities of baryons and radi¬ 
ation, respectively. The parameters Q b and Q r may be 
defined in the same way as in the ACDM model, namely 


n b = 


8nGp° 

~3Hl 


and fL = 


87 iGp® 

~3Hl 


(H10) 


where and are the densities at the present time 
and Ho is the present value of the Hubble parameter. 
Substitution into (H7) gives 


! 


y 3 + y 4 


^ + 2 a(HR/j\ 


3 V 




= K, 


(Hll) 


which has a unic^ue positive solution for 770 when a and b 
are non-negative. 

The integration was performed for different values of 
the parameters a, b and Q b and a good fit to the ob¬ 
served expansion history, given experimentally as a plot 
of residual magnitude against redshift for distant super¬ 
novae, was found for a range of values of a and b with 
fl b < 0.01. The time evolution of the scale factor in the 
string model for a = 0 , b = 0.01 and Q b = 0.01 is shown 
in Figure [i] The deceleration parameter, q = —RR/R 2 , 
is equal to —0.76 at the present time and the value of 
HR/j/c at the present time is 2.91. A similar time evo¬ 
lution can be obtained with Q b = 0.01 for other values 
of a and b , for example a = 0.15 and 6 = 0, for which 
q = —0.95 and HR/j/c = 2.79. The magnitude-redshift 
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FIG. 5: Plot of residual magnitude against redshift for the 
scale factor in Figure [2] with parameters a = 0, 6 = 0.01 and 
for the scale factor with parameters a = 0.15, b = 0, together 
with together with the experimental data points and error 
bars from the Supernova Cosmology Project and the High- 
z Supernova Search Team. Both sets of parameters give as 
good a fit to the data as the ACDM model. 


relations for both sets of parameters are shown in Fig¬ 
ure [5j together with the data from the Supernova Cos¬ 
mology Project 26] and the High-z Supernova Search 
Team [271 . 

The graph of logi? against logf in Figure [6] shows a 
transition from R ~ f 1 / 2 to R ~ t 2 / 3 at logy ss —2.4, 
corresponding to the transition from radiation domina¬ 
tion to matter domination at redshift z eq ss 250. 



FIG. 6: Plot of log R against log t for the best-fit scale factors 
R(t) in the string model. The scale factor is proportional to 

t l/2 

in the radiation era and t 2 ^ 3 in the matter era, exactly 
as in ACDM. 


It is remarkable that the evolution of the scale factor 
has the same form in both the radiation era and the mat¬ 
ter era as in the conventional ACDM model, for all values 
of a and b. The result may be understood by referring to 
equation (Hill. The quantity in square brackets varies 
very slowly compared with the factors of 1/y 3 and 1/y 4 
so, to a good approximation, R/j/R 3 is constant in the 
matter era and R/j/R 3 is constant in the radiation era. 


It follows from the definition (H8) that Ru is propor¬ 
tional to t when both R(t) ~ t 1/2 and when R(t) ~ f 2//3 , 
so R 2 v /R a is constant when R(t) ~ f 1 / 2 and R^/R 3 is 
constant when R(t) ~ t 2 / 3 . The string model therefore 
gives R(t) ~ t 1,/2 in the radiation era and R(t) ~ t 2 / 3 in 
the matter era, as claimed. 

It is also remarkable that, at least for some values of a 
and b , the string model explains the transition from decel¬ 
eration to acceleration. When t he tra nsition does occur, 
it is easy to show from equation (Hll) that the scale fac¬ 
tor has the form R ~ f 2 at late times when 6^0 and 
grows exponentially when 6 = 0. At late times, the ra¬ 
diation term in (Hll I is neglegible and HRjj/c increases 
with time, so it follows that R 3 /R 2 r is proportional to 
( HRjj/c ) 2 when 6 ^ 0 and proportional to HRu/c when 
6 = 0. Moreover, if R accelerates at late times then the 
integral in (H81 converges, so R.u has the same time de¬ 
pendence as R. Thus H ~ 1/y/R when 6^0, giving 
R ~ t 2 , and H ~const, when 6 = 0, giving exponential 
growth. 

The best-fit value of f \ in the string model is rather 
less than the value suggested by experiment, although 
the values suggested by different experiments are not yet 
in agreement. The smallest value of f \ consistent with 
nucleosynthesis calculations in the conventional model is 
fib = 0.015 [28] and the model in [25] with no dark matter 
gives fib = 0.016. 

If the value of H at the present time is taken to be 
70km/s/Mpc then the value of Rjj corresponding to a 
best-fit value of HRjj/c at the present time of 2.8 is about 
12Gpc, or about 3.6 x 10 26 m. Neglecting the contribu¬ 
tion from radiation, equations (H9) and (H10I give the 
relation 


GM V = ^(^)V 2 - (H17) 

Taking fib = 0.01 gives GMu = 0.04 Rjjc 2 and a corre¬ 
sponding value of Mu of about 10 22 Solar masses. 


Appendix I: ‘Dark matter’ in the string model 

Consider one of the Machian strings of a test mass 
m connected to a distant mass to. In the absence of 
any interactions between the Machian strings, the string 
is straight and has total energy Gmm/R , where R is 
the length of the string. The energy per unit length 
is T(R) = Gmm/R 2 , assuming the energy to be dis¬ 
tributed uniformly along the string. Now consider the 
interaction between the Machian strings of to and the 
Machian strings of a larger mass M. Let / be the change 
the energy per unit length of the Machian strings of to 
and suppose, as discussed in Section [V] that the value 
of / at a given point depends on the ratio of the den¬ 
sity of Machian strings of M at that point to the density 
of background strings. If A denotes the mass per unit 
length in a Machian string joining two particles of unit 
mass, the Machian strings within a distance r of the M 
contain a mass XrMMu■ The corresponding mass den¬ 
sity, Pm, in the Machian strings of M is therefore given 
by A-Kp M r 2 dr = XdrMMu, so that pm = XMMu/Anr 2 . 
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The mass density in the background Machian strings, 
Pb, is given by pb = {\M^Ru/2)/(A'kR^ j /3), so that 
PM/Pb ~ (M/M[/)(i?^/r 2 ). The ratio of densities is 
therefore given by the variable 


M /M[/\- 1 

|x - x M | 2 V Rfj / 


(II) 


The function / is assumed to be a function of u, so the 
energy per unit length of the Machian strings of m at the 
point x may be written in the form 


£(i?,x) = T(i?){l + /[n(x)]}. (12) 


At the position of the test mass m, the ratio of Machian 
string densities (Jll]) is equal to 


M / Mjj x- 1 
r 2 V R^j ) 



(13) 


where r is the distance between m and M, since ao ~ 
QAc 2 /Ru and, from Appendix[Hj GMu ss 0.047?[/c 2 . In¬ 
teraction between the Machian strings is therefore ex¬ 
pected to lead to corrections to the Newtonian law of 
gravitation with an acceleration scale of order O.lao- 
The magnitude of the expected additional acceleration 
can also be estimated by considering the limiting case of 
maximum asymmetry in the Machian strings of m, where 
all the strings point in the same direction at the centre. 
The tension T(R) = Grrim/R 2 in a string of length R 
connected to a distant particle of mass m makes a con¬ 
tribution to the total acceleration of magnitude Gm/R 2 . 
If all the Machian strings around m were arranged to 
point in the same direction at m, without changing their 
tensions, then the total acceleration would be given by 
integrating Gm/R 2 over all distant matter. The resulting 
acceleration would be 


/ G p 

dV = AitGpRu 
2 kc 2 


Ru 


since (B2) gives 


kc 2 = 


j ^dV = 2-KGpR'l,. 


(14) 


(15) 


In terms of the MOND acceleration scale qq ss 0.4c 2 /Ru, 
equation 


gives, with k = 3/40 from (B14|, 


0.4oq . 


(16) 


To calculate the additional acceleration explicitly, a 
particular form must be chosen for the function f(u). 
The limit u —> 0 corresponds to the limit of no interac¬ 
tion, so clearly f{u) —> 0 as u —>• 0. Appendix 12 shows 
that the condition / < 1 must be applied to ensure that 
the string tension remains positive, so f(u) must tend to 
some constant A < 1 in the limit u —> oo. To obtain flat 
galaxy rotation curves, the additional acceleration is re¬ 
quired to be proportional to 1/r at large distances from 


the centre of the galaxy. It turns out that a 1/r acceler¬ 
ation is obtained if f{u) ~ y/u as u —> 0, so consider the 
function fi(u) defined by 


h(u) 



(17) 


The corresponding Machian string paths were calculated 
using the method described in Appendix |I 1| The pattern 
of Machian strings for a test mass at the edge of a typical 
galaxy is shown in Figure [7J where the edge of the galaxy 
is defined as the distance, R g , at which the Newtonian 
acceleration is equal to ao- Since a o ~ 0Ac 2 /Ru, R g is 
given by 


GM 


0.4 


R u 


=> Rg 


0.32 R v 



(18) 



FIG. 7: Asymmetry of Machian strings for a test mass at 
the edge of a galaxy. The centre of the galaxy is represented 
by the black dot and the centre of the test mass is on the left, 
at the distance R g defined by equation (Jl8j) . All the strings 
shown are Machian strings of the test mass. 


In Appendix [I2j the tensions in the Machian strings at 
m are calculated and integrated over the asymmetrical 
Machian string distribution to find the total additional 
acceleration of the centre. In Appendix[[3j the additional 
force between the two masses is found independently by 
calculating the change in total energy of the Machian 
strings as a function of the distance between the masses 
and the two methods are found to be in agreement. For 
the function f 2 {u), defined by 


f2{u) 


Asju u < 1 
A u > 1, 


(19) 


which has the same limits as fi for both a < 1 and 
u 1, the formula for the acceleration calculated in 
Appendi x |I 3| can be evaluated analytically. 

Figure |8] shows the accelerations corresponding to the 
interaction functions /i and f 2 calculated in Appendix |T3| 
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FIG. 8: The additional accelerations, in units of ao, for the 
interactions between Machian strings corresponding to the 
functions /i and /> in equations © and ( [19] ), respectively, 
for A = 1. The green curve shows the average of the two 
accelerations. The radial distance r is in units of R g . 


in the limit A —> 1. For both curves, the additional ac¬ 
celeration for r R g is approximately a^Rg/r, which is 
exactly the asymptotic behaviour in MOND. The aver¬ 
age of the two additional acceleration curves, shown in 
green, is the curve used to plot Fig [3} 

It is important to check that the additional accelera¬ 
tion is negligible in the Solar System since no dark matter 
effects on the motion of planets or asteroids have been 
observed. Indeed, the total mass of any dark matter 
within the orbit of Saturn is estimated [3D] to be less than 
1.7 x 10 -10 Solar masses, or about 3 x 10 2O kg. The cor¬ 
responding Newtonian acceleration is 8.9 x 10~ 15 m/s 2 ss 
7.5 x 10 _5 a o , taking the orbit of Saturn to have radius 
r ~ 1.5 x 10 12 m, which provides an upper limit for the 
magnitude of the additional acceleration. The value of u 
at the radius of Saturn is about 4.4 x 10 5 , taking Mu = 
10 22 Solar masses and Ru = 3.6 x 10 26 m. According to 
the calculations in Appendices |I 2| and |I 3| below, the ad¬ 
ditional acceleration on a test mass m for the functions f± 
and /2 is approximately (3k/ y/uo)(c 2 / Ru) ~ 0.6ao/-y/uo 
for uq 3> 1, where Uq is the value of u at the centre of m. 
The corresponding additional acceleration at the radius 
of the orbit of Saturn is 9.0 x 10 -4 ao, which is too large 
by an order of magnitude, so it appears that a different 
function /(it) is needed. 

The functions /1 and / 2 defined above can easily be 
generalised by changing the value of u at which the tran¬ 
sition from / ~ Asjii for u <C 1 to / ~ A for u <C 1 oc¬ 
curs. Suppose, for example, that the length scale of /(it) 
is decreased by a factor A. The function / 2 (a) in ( [To] ) 
then becomes 


t , \ f Ay/Xu it < 1/A 

Mu) = \a 


( 110 ) 


The function / 3 (a) is larger by a factor of y/X for u <C 1, 
so the additional acceleration for Ho <C 1 is also expected 
to be larger by a factor of y/X. The effect on the addi¬ 
tional acceleration for ito 1 may be found by repeating 
the analytic calculation in Appendix |I3| W hen the func¬ 
tion / 3 is substituted into equation (1411, the function 


I(p) in equation (|I45|) becomes 


KH> 2 V 2 4 4p 


A 


' + y/X 


p - V A 

v P + V A > 


(Ill) 


where p = 1/yTto, which has the limits I(p) = \/A(l + 
In y f ^ — lnp + ...) for p > 1 and I(p) = \/A(l + 

In J -^/) — p 2 /(6\/A) + ... for p <C 1. The additional 

acceleration is therefore multiplied by y/X for p 1, 
i.e. for uq < 1, as expected, and divided by y/X for 
uq 1. Thus, for A > 1, the function / 3 (u) gives a 
larger additional acceleration in the MOND regime and a 
smaller additional acceleration in the Solar System. The 
exact form of the function f(u) remains to be derived by 
considering more carefully the nature of the interaction 
between the Machian strings. 


1. Calculation of Machian string paths 

Let s be the path length along one of the Machian 
strings of m connected to a distant mass to, with s = Si 
at to and s — Sf at to. If x(s) denotes the position of the 
point along the string at path length s then x(sj) = X; 
and x(s/) = X/, where X,; and X/ are the positions of 
the masses m and to, respectively. The total energy in 
the string is 


r s f 

E= £(R, x) ds , (112) 


where £ (R, x) is the energy per unit length defined in 
and R = Sf — Si is the length of the string. According 
to the IMC, the total energy of all the strings connected 
to the masses to and to is Es = me 2 + me 2 — E. To 
minimise the total energy of the system it is therefore 
necessary to find a string path for which the energy (1121 


is a maximum. Since the energy (112) tends to zero in 


the limit that the string becomes infinitely long, string 
paths that maximise (112) do exist. 


Consider a variation <5x of the string path with the 
string held fixed at the distant mass to, so that Sx = <5Xj, 
say, at s = Si and 6x = 0 at s = Sf. Since the total 
length of the string changes it is convenient to introduce 
the parameter a along the string path so that a is fixed 
at to and to, with a = 0 at to and cr = 1 at to. The 
energy E is then given by 


E = f £(R, x) |x'| da , 
Jo 


(113) 


where the prime denotes differentiation with respect to 
cr, and the string length R is given by 


R = 


da. 


(114) 
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The variation of (1131 is 


6e =1\(m sr + V£ ^) |x,| + f l^f} da ’ (I15) 

(116) 


where 


5R = d*. 

Jo l x l 


Integrating (116) by parts gives 
r 1 


/x'Axn 1 f f x" x'(x'.x")'| 

“ = (Wo“l 0 ( ’ 

After changing back to the path length parameterisation, 
for which |x'| = 1 and x'.x" = 0, (117) becomes 


SR = —x'(si).SX; — / x" .Sx ds. 


(118) 


Similarly, (115) becomes 


5E = — £(!?, Xi) x!(si).5Xi 


V£ — x'(x'.V£) — £x" Axj ds. (119) 


The variation SR is independent of s and may be taken 
outside the integral. Substituting for SR from (118) then 
gives 

5E = — |£( J R,X J ) + /(i?)]x , (s i )AX i 

+ J * -x'ix'.V£)-[£ + I{R)]x"YSxds, ( 120 ) 

r Sf Q£ 

where I(R) = / —— ds. (121) 

J Si ('JA, 

The requirement that E is stationary for all variations 
of the string path for which the string is fixed at both 
ends, i.e. for which 5X; = 0, gives the path equation 


V£ - x'(x'.V£) -[£ + J(i?.)]x" = 0. 


( 122 ) 


After substituting x" = m, where k is the curvature 
and n is a unit vector normal to the tangent vector x!, 
the component of equation (122) along x' is found to 


be identically zero and the component along n gives the 
equation for the curvature, 


c(x) = 


n.V£ 


£(R,x) + I(R) ' 


(123) 


A computer program was written to calculate the paths 
of Machian strings around a test mass to, with the energy 
per unit length of the Machian strings defined by ( fl2j ) and 
the interaction function /i(u) defined by |l7|. If ip de¬ 
notes the angle that a Machian string makes with the 
line of centres joining to and M, the path equations are 
x' = cost/), y' = sin ip and ip' = k, where the prime de¬ 
notes differentiation with respect to path length and n is 
given by equation (123). The value of A was taken to be 


0.5 to avoid possible numerical problems associated with 
vanishing string tension when A is close to 1. Machian 
string paths were calculated for different values of A/, for 
a given separation, r, of the two masses. Figure [T] shows 


a two-dimensional set of string paths, lying in a plane 
containing the line of centres, for the case r = R g . Since 
the effect of the mass M is negligible at large distances, 
the initial directions of the strings were adjusted to en¬ 
sure that the distribution of strings at large distances is 
uniform. The two-dimensional set of string paths then 
represents a cross-section of the actual three-dimensional 
Machian string distribution. The Machian strings around 
to are seen to be asymmetric, with a higher density of 
strings on the side nearest to M. 


2. Direct calculation of the additional acceleration 


For paths satisfying the path equation (122), it follows 


from (120) that the change in the total energy of the 


string when mass to is displaced by <5X; is 


SE = - 


£(R,Xi) + I(R) x , (s i ).6X i . 


(124) 


The total energy of the system is Es = me 2 + me 2 — E, 
so 


SE s = 


£(R,Xi) + I(R) x'(si).SXi. 


(125) 


If F denotes the force exerted by the string on the mass 
to then the work done by the system is FAX;, so SEs = 
—FAX,-. The force exerted on the mass to is therefore 


F = 


£{R,X t ) + I{R) x'(si ), 


(126) 


where x' is the unit vector from to to to. The string 
tension at a general point x along the string is given by 


T(R,x) = ~[£(R,x)+I(R)] . 


(127) 


Note that, in the absence of any inter actio ns between the 
strings, £ = T(R) = Gmm/R 2 and (126) reduces to the 
Newtonian gravitational force 


F = 


Gtoto 


(128) 


Consider the energy per unit length £ (R, x) d efined 
by equation p2]). After substituting ([12]) into (121), not¬ 
ing t hat T(R) is proportional to 1/Fir, the string ten¬ 
sion (127) becomes 


T(R,x) = T(R)F{R,x), 


(129) 


say, where 


F(R,x) = 1 - f[u{x)] + — / f[u(x(s))]ds. (130) 


To ensure that the curvature (123) remains finite, the 


string tension must be positive everywhere along the 
string so the function F(R , [x]) must be positive. For 
a stri ng of le ngth Rjj, the change of variables s = 
a Ru y/M/Mjj gives 


F(R,x) = l-/[u(x)]+2 A /— ^ 


M 


f(u ) da , 


( 131 ) 
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where u = (p 2 + a 2 — 2pcrcos0) _1 , p is defined by 
r = pR\j yj M/Mjj and 9 is the angle between x and 
xjvf. For the functions f(u ) defined in © and p9| , 
/(zt) ~ y/u ~ 1/cr for cr 1. The contribution to the 


integral in (131) is of order unity from a < 1 and of order 


In (Mjj/M) from a > 1, so the third term in ( |I31[ ) is of or¬ 
der yj M/Mir In (Mjj/M) and is therefore negligible since 
Mjj/M is very large. It follows that F(R,x) < 1, so the 
interaction with M reduces the tension in the Machian 
strings of m. The condition needed to ensure that the 
string tension remains positive is / < 1. 

For a given two-dimensional set of Machian string 
paths, such as those shown in Figure [7J the additional 
acceleration acting on the test mass to may be calcu¬ 
lated as follows. Since the distribution of strings at large 
distances is uniform, the total force acting on the centre 
is given by 

F = y J T ( 0 ) sin0d0 # = \ j T(0)sin 0d6, 

(132) 

where T(0) is the tension at the origin in a string whose 
direction at large distances is at an angle 9 to the line of 
centres. When A < 1, the integral in (131) is negligible 


and the magnitude of the string tension at m is T(R)[1 — 
/(zio)] for all strings, where uq is the value of u at to. The 
component of the total force along the line of centres, 
from M to m, is then 


F = —[1 — /(«o)] J cosipdcos9, 


(133) 


where ip denotes the angle made by the initial direction 
of a string, at the centre of to, with the line of centres. 
According to equation ([l6]), the acceleration of to due 
to the Machian strings is equal to 2kc 2 /Rjj when the 
tension in all the strings is T(R) and all the strings point 
from to to M at the centre of to, i.e. when ip = 0 for 
all strings. The a ccele ration due to the Machian strings 
corresponding to (133) is therefore 


9kt 2 

a = —[1 — f(uo)](cosip ), 


Ru 


(134) 


where (cos ip) is the average value of cosz p, as a function 
of cos 9 , when cos 9 uniformly distributed in the interval 
[- 1 . 1 ]- 

The quantity (cos ip) defines the asymmetry of the dis¬ 
tribution of Machian strings. The asymmetry is zero for 
M = 0, when the Machian strings of to are spherically 
symmetric, and unity in the limit that all strings of m 
initially point towards M. A plot of cos ip as a function 
of cos 9 for the set of Machian string paths in Figure [7] 
is shown in Figure [9] The distribution of strings and 
the resulting asymmetry was calculated for different val¬ 
ues of M, corresponding to values of uq in the range 
10 -6 < uq < 10 6 . It may be seen from Figure 10 that 
the asymmetry is largest for zto ~ 1. The variation of the 
asymmetry for zto <C 1 and zto 1 is given by 


(cos ip) 


f 0.5^u 0 


u, 0 for zt 0 <C 1 

Z/y/uQ for u 0 » 1. 


(135) 



FIG. 9: A plot of cos ip as a function of cos# for the strings 
in Figure[7j where ip is the angle that a string makes with the 
horizontal and 9 is the angle of the string at a large distance 
from M. The asymmetry, given by averaging the function 
over the range [—1,1], is equal to 0.16. 



FIG. 10: A log-log plot of the asymmetry in the Machian 
strings of a test mass due to the presence of a mass M at a 
distance r as a function of zto, where zto is the ratio of M/r 2 
to Mu/Ru- 


The corresponding acceleration is cal culat ed for the data 
points in Figure 10 using equation (134), and is shown 
in Figure [IT] The additional acceleration for the function 
/i(zt), with A = 0.5, therefore has the limits 


2 kc 2 j 0.5^/zip for zt 0 1 
Ru \ 1-5/v^o for zt 0 »l. 


(136) 


For comparison with the results of Appendix |I 3[ it is 
convenient to rewrite the limits in terms of the variable 
p defined by r = pR\j\/M/Mu , which is related t o zto 
by the simple equation uq = 1/p 2 . In terms of p, (136) 
becomes 


2 kc 2 ( 1.5 p for pCl 
Ru \ 0.5 /p for p » 1. 


(137) 
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Changing variables from cos 8 to u , using du = 
2 pau 2 d cos 6 , gives 

~ kwc 2 MRu f V da /(zt) 

~ 2r M U Jo f 7 J(p+C t)-2 M 2 

(141) 

When the IMC is taken into account, the corresponding 
change in energy of the whole system, including all the 
other masses other than to and M, is equal to —A E. 
The additional acceleration of the mass m due to the in¬ 
teraction between the Machian strings is therefore given 

by 


FIG. 11: The magnitude of the additional acceleration, in 
units of c 2 / Ru, calculated from the asymmetry of string paths 
(green crosses) using the data points in Figure |10[ for the 
interaction defined by the function (|T7[ ) with A = 0.5. The 
blue curve is calculated in Appendix |I 3| using an approximate 
formula for the total interaction energy. 


3. Calculation of the additional acceleration from 
the Machian string interaction energy 


The change in the total energy of the Machian strings 
of m corresponding to the interaction with the Machian 
strings of M defined by equation © is 

ae = Gl ft™ 1 j o /[ u ( x )W s > ( I38 ) 


where s denotes the path length along the i th string. Al¬ 
though the Machian strings around to are asymmetrical, 
as illustrated in Figure [TJ a symmetrical distribution of 
strings around m may nevertheless be assumed in order 
to calculate an approximate formula for AE. The total 
mass to in an elemental solid angle dfl and thickness dR 
at radius R is pR 2 dRd£l 1 where p is the average matter 
density, so (138) then becomes, after integrating over the 
azimuthal angle, 


pRu /‘ 7r pR 

AE S3 27 xGmp / dR sin d dd / /[w(x)]ds. (139) 


/ o 


lo 


/ 0 


Let r be the distance between the centres of the two 
masses, so that |x — xm| 2 = r 2 + s 2 — 2r scos6*. Since 
/ —> 0 as s —> oo, the integral over s in (139) is insensitive 


to the value of R, so R can be replaced by Ru in the 
upper limit of the integral over s and the integral over 
R then simply gives a factor of Ru- If the dimensionless 
variables p and a are defined by r = pRu \ JM/Mu and 
s = aRu \/MJ Mu then u = (p 2 + a 2 — 2pa cosf?)” 1 . 
Equation (139) then becomes, using 


AE 




f[u(x)] dcosd. 


1 dAE 
to dr 


(142) 


For the function fi(u), the energy (141) can be eval¬ 


uated numerically and the result, for A = 0.5, is shown 
as the blue curve in Figure [ll] The corresponding limits 


are 


2 kc 2 ( 1.5 p p < 1 

\ 0-4 /p p > i, 


(143) 


which is seen to be in good agreement with (137) and 
confirms the validity of the above calculations and the 
calculations in Appendix |l2] with an accuracy of ~10%. 



FIG. 12: The magnitude of the additional acceleration, in 
units of c 2 /Ru, for the interactions between Machian strings 
corresponding to the functions fi and in equations © 

and fl9| ), respectively, for A = 0.^ _ Both calculations are 

based on the approximate formula (1391 for the total interac¬ 
tion energy. The blue curve is the same as in Figure © 


If f(u) is take n to be the function / 2 (u) defined in 
the energy (141) can be evaluated analytically. The result 


is 


AE ss 2 /cAtoc 2 


M 

Mu 


I(P), 


(144) 


where the function I(p) is given by 

p+i , ln Mu/m ^' 


I{p) = h(l- E 4-u} <u 


P~ 1 


( 145 ) 


( 140 ) 
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and the corresponding additional acceleration is 


a ~ 2nA—— I' (p) ■ 
Ru 


(146) 


For p » 1, I(p) = 1 + In \ — lnp + ... and the 


additional acceleration is therefore A/p times c 2 /Ru , 
which is indeed proportional to 1/r. For p <C 1, I(p) = 


1 + In J — p 2 /6 + ... so, for A = 0.5, (146) has the 


limits 


2kc 2 ( p/6 for p < Cl 
Ru \ 0.5 /p for p > 1. 


( 147 ) 


Waves with angular frequency to and amplitude A have 
energy F A 2 oj 2 /2£ 2 c 2 per unit length along the string, 
so the energy flux in the string is F A 2 uj 2 /2t/ 2 c g . If all 
the Machian strings of the mass m, : are assumed to have 
waves of amplitude Ai, the total rate of gravitational ra¬ 
diation from rrii is 


E 


GmimjAfoj 2 

2r ij c g 


SGrriiMuA/tu 2 

2 C gRfj 


(J3) 


With ui 2 = GM/r 3 , the total rate of energy loss due to 
gravitational radiation in the strings connected to the 
two masses is therefore 


The accelerations corresponding to the two functions fi 
and f 2 are shown in Figure [l2| 

To plot the curves shown in Figure [8j the values on 
the y-axis were multiplied by a factor of 2.5 to convert 
from units of c 2 /Ru to units of a 0 . The values on the 
x-axis were multiplied by a factor of 3.2 to convert from 
p to the distance r measured in units of R g since, from 
equation |k]), r/R g = 3.2 p. 


Appendix J: Gravitational radiation 


Consider a binary pulsar consisting of two masses, mi 
and m 2 , in orbit around each other. For the special case 
when the masses move in a circular orbits, the masses 
have constant angular velocity ui = \JGM/r 3 , where 
M = 777 1 + 7772 and r is the separation of the two masses. 

In the string model, each mass exerts a periodic dis¬ 
tortion in the Machian strings of the other as it or¬ 
bits around it. The distortions have angular frequency 
ui and are expected to generate gravitational waves of 
angular frequency w in the Machian strings of each 
mass. The speed of the waves can be calculated from 
the Lagrangian ( |A6| ) by considering perturbations in a 
Machian string of the form X 11 = (ct,£<j + <5x). Then 
{dX^/dT)(dX^/da) = -£.<5x - «5xAx', (dX”/da) 2 = 
—t/ 2 — 2£.(5x' — (Sx' 2 and s = 1 — <5x 2 /c 2 . For trans¬ 
verse perturbations, with £.<5x = 0, expansion of (A6) to 
quadratic order gives 


cLE __ 3G 2 M u {m 1 Al + m 2 A\)(m.i + m 2 ) , , 

dt 2 c g R 2 r r 3 


The amplitudes A\ and A 2 of the gravitational waves 
in the strings of the two masses could be found by gener¬ 
alising the numerical calculation of the field string paths 
in Appendix ||J The interaction between the Machian 
strings of mi and m 2 has so far been considered only 
for the case when mi is very much smaller than m 2 , so 
the interaction must first be defined for the case when the 
masses are of comparable magnitude. The paths of the 
strings would also have to be calculated as a function of 
time and such calculations have not yet been attempted. 

The prediction of General Relativity |3I] > as confirmed 
experimentally by observations of the decay of the orbital 
period of binary pulsars [321 EH], is 


dE 32G 4 (mim 2 ) 2 + m 2 ) 

dt 5r 5 c 5 


Comparison with equation (J4) shows that the amplitude 
of the waves in the Machian strings of 7Tii required for 
agreement with experiment has the form 



with a similar expression for A 2 , where the dimensionless 
parameter 77 is given by 


c = r (i + oi^ + ...)(f 

F / Sx 2 dx' 2 


2 


+ 5x' 2 + ... 


- 1/2 


(Jl) 


The equation of motion for <5x corresponding to © 
is 2ai(5x/c 2 = Sx .”/£ 2 , so that transverse gravitational 
waves propagate along the strings with speed 


c 


(J2) 


r/ 2 = 


32 

15\/2a/ ‘ 


(J 7 ) 


The value of 77 corresponding to the value of ai given 
by (B13) in Appendix 11^ is 77 = 0.584. In equation (J 6 ), 


the length scale Ru\Jm\/Mu is the distance around 7771 
within which the density of Machian strings of 7771 is 
larger than the density of background strings and the di¬ 
mensionless factor Gm 2 /rc 2 is a measure of the strength 
of the gravitational interaction at mi due to m 2 . 
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